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' We use the microscopic Hartree-Fock approximation to investigate various quantum phase tran- 



sitions associated with possible spontaneous symmetry breaking induced by a tilted magnetic field 
in the integral quantum Hall regime of wide parabolic wells and zero width double well (bilayer) 
systems. Spin, isospin (associated with the layer index in the bilayer systems), and orbital dynamics 
all play important roles in the quantum phase transitions being studied. We propose a general class 
of variational wavefunctions that describe several types of parity, spin, and translational symmetry 
breaking, including spin and charge density wave phases. In wide well systems at odd filling factors, 
we find a many-body state of broken parity symmetry for weak in-plane magnetic fields and an 
isospin skyrmion stripe phase, which simultaneously has isospin and charge modulation, for strong 
in-plane fields. In wide well systems at even filling factors, we find direct first order transitions be- 
^ , tween simple (un)polarized QH states with several complex many-body states that are only slightly 

higher in energy (within the Hartree-Fock theory) than the ground state. We suggest that going 
beyond the approximations used in this paper one may be able to stabilize such many-body phases 
^ . with broken symmetries (most likely the skyrmion stripe phase). In a bilayer system at the filling 



factor V — 4:N ± 1, where A'' is an integer, we obtain an isospin spiral stripe phase in addition 
to the known (in) commensurate phases and the stripe phase without isospin winding. We do not 
find a charge or spin density wave instability in the bilayer system a,t u = AN + 2, except for the 
' known commensurate canted antiferromagnetic phase. Zero temperature quantum phase diagrams 

f— I , for these systems are calculated in the parameter regime of experimental interest. We discuss the 

Q ' symmetry properties of our predicted quantum phase diagrams and give a unified picture of these 

O ' novel many-body phases. We point out how quantum level crossing phenomena in many situations 

(tuned by the applied tilted magnetic field) may lead to interesting quantum phases and transitions 
among them. A conceptually new aspect of our theory is the predicted possibility for the spon- 
^ ] taneous breaking of parity symmetry, which indicates a "ferroelectric" quantum order in integer 

. quantum Hall systems and has not been considered in the literature before. 
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O . I. INTRODUCTION 

Unidirectional charge density wave order (also called stripe order) in quantum Hall (QH) systems has been ex- 
tensively studied [1] since the first theoretical prediction in 1996 [2,3] and the first experimental observation in high 
Landau levels via the magnetoresistance anisotropy measurement in 1999 [4,5]. Many related phenomena, e.g. trans- 
^ port via internal edge state excitations [6-8], liquid crystal phases [1,9], reorientation of stripe directions [10-12], 
O , and re-entrant integer quantum Hall effect [13], have been widely explored both theoretically and experimentally in 
O ■ this context. Stripe formation at fractional filling factors, v = (27V -|- 1)/(47V -|- 4), corresponding to the composite 
^ , fermion filling factor, v* — N + 1/2, was also proposed to have energy lower than the conventional Laughlin liquid, 
composite-fermion Fermi sea or paired composite- fermion state [14]. (Throughout this paper TV is zero or an integer, 
= 0, 1 • • •, and = 1, 2, 3 • • • is the Landau level filling factor of the whole system.) However, most of the stripe 
^ . phases discussed in the literature so far are formed by electrons in the top (half-filled) level of the QH system at high 
■ ■ ' half-odd-integer (HOI) filling factors (e.g. i/ = N +1/2 for N > 4), while the stripe formation at integer filling factors 
is seldom studied either theoretically or experimentally. 

We develop a detailed theory in this paper for possible spontaneous symmetry breaking and associated exotic 
quantum order in wide-well and double-well integer quantum Hall systems by considering all the symmetry proper- 
ties of the realistic system Hamiltonians. The great advantage of the quantum Hall systems in studying quantum 
critical phenomena is, however, the existence of various energy gaps (at the Fermi level) which enables us to carry 
out reasonable energetic calculations (within a mean-field Hartree-Fock theory) to quantitatively check whether the 
various possible quantum phase transitions and exotic quantum order allowed by symmetry considerations are actu- 
ally energetically favored in realistic systems. We therefore construct explicit variational wavefunctions for various 
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possible (exotic) quantum states, and carry out energetic calculations to find the optimal ground state. Using an 
in-plane parallel magnetic field (in addition to the perpendicular field necessary in the quantum Hall problem) to tune 
system parameters, we find a surprisingly rich quantum phase diagram in our systems of interest. We believe that 
the symmetry-broken states with exotic quantum order predicted by our theory should be experimentally observable 
in transport and inelastic light scattering experiments. We emphasize that the salient feature of our work (mak- 
ing it particularly significant from the experimental perspective) is that we not only identify spontaneous symmetry 
breaking and exotic quantum order allowed by the system Hamiltonian, but also carry out Hartree-Fock energetic 
calculations to obtain the quantum phase diagrams using realistic system Hamiltonians. Among the spontaneous 
symmetry breakings predicted in this paper arc, in additional to the usual translational and spin symmetries, the 
spontaneous breaking of discrete parity symmetry in a number of interesting situations. 

Very recently, magnetoresistance anisotropy was observed in both doped GaAs/AlGaAs [15] and Si/SiGe [16] based 
semiconductor quantum well systems at even integer filling factors, when a strong in-plane magnetic field (B||), in 
addition to the perpendicular magnetic field {B±) producing the 2D Landau levels, is applied. In Refs. [15,16], it is 
proposed that the strong in-plane magnetic field increases the electron Zeeman energy so much that the highest filled 
Landau level (with spin down) has a level crossing with the lowest empty Landau level (with spin up), and then a 
charge density wave (CDW) phase may be stabilized by electron Coulomb interaction. Our recent work [17] on the 
magnetoplasmon energy dispersion of a wide well system in the presence of a strong in-plane magnetic field finds a 
near mode softening in the spin-flip channel, suggesting the possibility of a spin density wave (SDW) instability near 
the degeneracy point. As a result, we proposed a spin-charge texture (skyrmion [18]) stripe phase [8,17] to explain the 
magnetoresistance anisotropy observed in the experiments [15,16]. In an earlier theoretical work, Brey [19] also found 
a charge density wave instability in a wide well &i v = 1 without any in-plane magnetic field. Murthy [20] recently 
considered the the coexistence of quantum Hall order and density wave order in a single well at = 2 with zero 
Zeeman energy and strong level mixing. To the best of our knowledge, however, systematic theoretical analyses of 
these candidate stripe phases at integer QH systems have not yet been carried out and need to be further developed. 

A double well (bilayer) quantum Hall system is another system where one may see a stripe phase order at integral 
filling factors. At total filling factor v = AN ± 1, electrons are equally distributed in the two layers if no gate voltage 
or in-plane magnetic field is applied. We note that the bilayer systems of filling factor u = AN — 1 {N > 0) can 
be understood to be equivalent to v = 4{N — 1) -|- 1 = 4 A'' — 3 system by electron- hole symmetry in the top filled 
Landau levels, if the Landau level mixing is negligible. Therefore all of our results shown in this paper for bilayer 
1/ = AN + 1 quantum Hall systems can be applied to v = A{Nj^ 1) — 1 = 4A'' — 3 system as well. We will not distinguish 
these two systems and will mention v = AN + 1 case only in the rest of this paper. When the layer separation, rf, 
is small compared to the 2D magnetic length, spontaneous interlayer coherence can be generated by interlayer 
Coulomb interaction even in the absence of any interlayer tunneling [21,22], and the ground state is therefore best 
understood as a Halperin (1,1,1) coherent phase with finite charge gap [23]. When d is comparable to or larger 
than however, there should be a transition to two decoupled compressible v = 2N + 1/2 systems or one of the 
competing many-body phases. A detailed analysis of all competing phases is still lacking, although it is commonly 
believed that for = there is a direct first order transition from the (1,1,1) phase to the compressible states, 
although for A'' > there may be intermediate quantum Hall phases with stripe order [24] . Other more exotic phases, 
including the ones with electron pairing [25], have also been proposed in the literature. This interlayer coherence 
and the corresponding Goldstone mode of the symmetry-broken phases have been extensively studied theoretically 
[26-28] as well as experimentally [29,30] in the literature. In addition to the layer separation, the in-plane magnetic 
field is another important controlling parameters for the double well QH systems. In the presence of an in-plane 
magnetic field, an Aharonov-Bohm phase factor, associated with the gauge invariance, has to be considered for the 
electron tunneling amplitude between the two layers and may cause a commensurate-to-incommensurate (C-I) phase 
transition at a critical magnetic field strength [26,31,32]. Such C-I phase transition arises as a result of competition 
of the interlayer tunneling and Coulomb interlayer exchange interaction. The former favors a commensurate phase, 
in which the phase of the order parameter winds at a rate fixed by the parallel magnetic field; whereas the latter 
favors an incommensurate phase, in which the order parameter is more nearly uniform in space [22]. Taking into 
account the stripe order associated with the layer separation (d) and the spiral order associated with the in-plane 
magnetic field (By) the possibility exists for a very complex and rich phase diagram with many competing phases 
occurring for various values of d and in the bilayer double well system a,t v = AN + 1. Following our earlier work 
based on the collective mode dispersion [8,17], in this paper we carry out the ground state energetic calculation within 
the Hartrce-Fock approximation to obtain and describe these exotic phases, which break isospin rotation, parity, 
and/or translational symmetries. We will also investigate the possibility of stripe formation at the even filling factor, 
1/ = AN + 2, where a nonstripe canted antiferromagnetic phase (CAE) has been proposed and extensively studied in 
the literature [33 41]. The goal of our work is to discuss possible symmetry broken exotic quantum phases in bilayer 
and wide- well quantum Hall systems. 

In this paper, we use the idea of isospin to label the two (nearly) degenerate energy levels at the Fermi energy, 
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and construct trial many-body wavcfunctions that have various kinds of (possibly exotic) isospin order. The relevant 
isospin quantum number can be associated with Landau level index, layer index, spin index, or any other good quantum 
numbers describing the corresponding noninteracting system. Within Hartree-Fock approximation we investigate 
phases with isospin spiral and/or stripe orders and discuss their relevance for various systems. In Table I we summarize 
the definition of the isospin up(down) components for each system individually. For the convenience of discussion, 
we use following notations and labels to denote the systems in the rest of this paper (see also Fig. 2): VKl(Ty2) and 
VF1'(VF2') denote the intersubband and intrasubband level crossing (or level near degeneracy in Wl' case) in a wide 
well system with filling factor v = 2N + 1{2N + 2) respectively; Dl and D2 are for the level crossing (or level near 
degeneracy in Dl case) in a double well system at total filling factor, jy = 4N + 1 and 47V + 2, respectively. (Note 
that the distinction between intrasubband and intersubband level crossing of a wide well system in the presence of 
an in-plane magnetic field is somewhat ambiguous and will be clarified in the next section.) By "level crossing" we 
mean a point in the parameter space of noninteracting electrons, for which the energy difference between the highest 
filled Landau level and the lowest empty Landau level vanishes, while by "level near degeneracy" we mean a region in 
the parameter space, in which the energy difference between the above noninteracting electron levels is not zero, but 
relatively small compared to the interlevel interaction energy, which can strongly mix the two close noninteracting 
Landau levels in certain situations. In this paper, we consider the following six quantum phases generated by our 
trial wavefunction: conventional incompressible integer QH liquid (i.e. isospin polarized), isospin coherent, isospin 
spiral, isospin coherent stripe, isospin spiral stripe and isospin skyrmion stripe phases, classified by the behaviors of 
the expectation values of their isospin components, {I±) and (J^). More precise definitions and the related physical 
properties of these phases will be discussed in Section III. In the last two columns of Table I we list the exotic 
many-body phases that become the ground states near the appropriate degeneracy points and the associated broken 
symmetries. For a wide well system a,t v = 2N + 1 near the intersubband level crossing (i.e. Wl in small in-plane 
field region), we find that interactions stabilize the isospin coherent phase, which breaks the parity (space inversion) 
symmetry and can be understood as a quantum Hall ferroelectric state [42] . For the same system near the intrasubband 
level near degeneracy region (i.e. Wl' in large in-plane field region), we obtain the isospin skyrmion stripe phase, 
which breaks both parity and translational symmetries of the system. For wide well systems at even filling factors 
we find that many-body states with isospin coherence and/or stripe order are not favored at either the intersubband 
{W2) or the intrasubband {W2') degenerate points within the HF approximation, exhibiting instead a trivial first 
order phase transition at the level crossing points. However, we find the calculated HF energy difference between 
the uniform ground state and the phases of spiral and stripe orders to be rather small, indicating the possibility of a 
spontaneous symmetry breaking in more refined approximations. Therefore, based on our HF calculation results, we 
suggest the existence of skyrmion stripe phase, breaking both spin rotational and translational symmetries near the 
degenerate point of W2' systems, may be responsible for the resistance anisotropy observed recently in Ref. [15]. (As 
for the magnetoresistance anisotropy observed in Si/SiGe semiconductor [16], the complication of valley degeneracy 
in Si makes a straightforward comparison between our theory and the experiments difficult.) For a double well system 
at z/ = 4A'' -|- 1 (Dl), several interesting phases can be stabilized in the parameter range of interest (see Table I), 
while only a commensurate CAF phase is stabilized for even filling factor, i/ = AN + 2 {D2), being consistent with 
the earlier results [36]. All of our calculation are carried out at zero temperature without any disorder or impurity 
scattering effects. We suggest that these symmetry broken phases discussed in this paper should be observable in 
transport and inelastic light scattering experiments. 

We note that the presence of an in-plane magnetic field -B|| (in the 2D plane) automatically introduced a tilted 

magnetic field of amplitude Btot ~ y^-Bjj -|- Bj^ at an angle of tan~^(B||/S_L) tilted with respect to the normal vector 

of the 2D plane since the field B±, perpendicular to the 2D plane, is always present in order to produce the 2D 
Landau level system. The presence of B|| (or more generally, the tilted field) introduces qualitatively new physics to 
2D quantum Hall systems by allowing a tuning parameter (i.e. _B||) without affecting the basic Landau level structure. 
Changing in a continuous manner may enable the 2D system to undergo various quantum phase transition which 
would not otherwise be possible. Thus the physics of quantum Hall systems becomes considerably richer in the 
presence of a continuously tunable parallel (or tilted) magnetic field. The general (and somewhat ambitious) goal of 
our theoretical work present in this paper is to study this rich quantum phase diagram of 2D quantum Hall systems 
as a function of the in-plane field (used as a tunning parameter). As such we concentrate in the two most promising 
candidate systems, namely the (single) wide quantum well system and the (bilayer) double quantum well system, where 
the tilted field configuration is expected to give rise to exotic symmetry-broken quantum orders in the quantum Hall 
ground states. The underlying idea is to use the in-plane magnetic field as a tuning parameter to cause level crossing 
(or almost level degeneracy) in the noninteracting system (around the Fermi level) to see if the interaction effects then 
drive the system to non-trivial symmetry-broken many-body ground states (particularly with exotic quantum order) 
which can be experimentally studied. Earlier examples of such exotic symmetry-broken ground states (arising from 
level crossing phenomena) include, in addition to our prediction [17] of a skyrmion stripe phase, the bilayer canted 
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antiferromagnetic phase predicted in Ref. [33] where a tilted field is not necessary. 

This paper is organized as follows: in Section II we introduce the Hamiltonians of interest, for both wide well and 
double well systems in the presence of a finite in-plane magnetic field. In Section III we propose various trial many- 
body wavefunctions incorporating both isospin stripe and isospin spiral orders. Physical features of the six typical 
phases generated by the wavefunctions will also be discussed in detail. For the systems of even filling factors, we 
propose trial wavefunctions that involve the four Landau levels closest to the Fermi energy simiiltaneously. (Including 
four rather than two levels in many-body wavefunctions was shown to be crucial for establishing the many-body canted 
antiferromagnetic state in bilayer systems at u = 2 [33].) We will then divide the Hartree-Fock variational energies 
and the related numerical results into the following four sections: wide well systems aX v ~ 2N + 1 (Section IV) and 
at V = 2N + 2 (Section V); double well systems at total filling factor v = AN + 1 (Section VI) and at u = 4iV 4- 2 
(Section VII). For the sake of convenience (and relative independence of different sections), results of each section are 
discussed independently and then compared with each other in Section VIII, where we also make connections to the 
earlier works in the literature. Finally we summarize our paper in Section IX, and discuss some open questions. 

II. HAMILTONIANS 

In this section, we present the Hamiltonians of the systems we will study in this paper, including both the single 
wide (parabolic) well systems and the double (thin) well (or bilayer) systems. Most of the formulae given in this 
section exist in the literature (perhaps scattered over many publications) and hence we will not derive them in details 
unless absolutely necessary. To make the notations consistent throughout this paper, we use the superscript W to 
denote physical variables or quantities for the wide well system and the superscript D for the double well systems. 
We will also describe the symmetry properties of these systems, which are crucial in formulating the many-body 
trial wavefunctions in the next section and in discussing the nature of the spontaneous symmetry breaking and the 
associated exotic quantum order in our various proposed phases. 

A. Wide Well System 

For a wide well system we consider a parabolic confinement potential [17] in the growth direction [z), Uc{z) = 
^m*u)QZ^, where m* is electron effective mass and wq is the confinement energy. The advantage of using a parabolic 
well model is that we can easily diagonalize the noninteracting Hamiltonian even in the presence of in-plane magnetic 
field. We will use eigenstates of the non-interacting Hamiltonian as the basis functions for writing variational Slater 
determinant states for the many-body wavefunctions that may be stablized by Coulomb interactions. We believe 
that our wide well results obtained in this paper for parabolic confinement should be qualitatively valid even for 
nonparabolic quantum wells. To incorporate both the in-plane magnetic field (i?||), which we take to be along the x 
axis in this section, and the perpendicidar (along the z axis) magnetic field {B±) in the Hamiltonian, we choose two 
kinds of Landau gauges for the vector potential: in one of them, yl[j^](r ) = {0,B±x — i3||Z,0), particle momentum 

is conserved along y direction, and in the other one, A[^-^{r) = [—B±y,—B\\z,0), particle momentum is conserved 
along X direction. The final physical results are of course independent of the choice of the gauge. However, as will 
be clear from the discussion below, various phases may be more conveniently discussed in different gauges, since the 
physical or the mathematical description of particular many-body states may be more natural in specific gquges. For 
the consistency of notations, all the explicit calculations presented in this paper will be done in the gauge A[y-^{r). 
Since the generalization of these equations to the other case is straightforward (see Appendix A) , we will only present 
the final results for the A^^j (r ) gauge. The noninteracting (single electron) Hamiltonian in the parabolic potential 
with both perpendicular and in-plane magnetic fields in gauge A[y]{f) is (we choose U = l throughout this paper) 

i?o^ = ^ + — fp. + ^ - + — + --*-o-^ - (1) 

° 2to* 2rn* V c c / 2m* 2 ° ^' ^ ' 

where Sz is the z component of spin operator and uj^ is Zeeman energy, proportional to the total magnetic field, 

Btot = \J + It been shown [17] that the one electron energies and wavefunctions of Eq. (1) can be 

obtained analytically even in the presence of an in-plane magnetic field by rotating to a proper coordinate. The 
resulting noninteracting Hamiltonian is a sum of two decoupled one-dimensional simple harmonic oscillators with 
energies [17], 
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Wl 2 = 



(2) 



i^n + wif is the effective 



where oj±^\\ = eB±^\\/m*c (i.e. is the conventional cyclotron frequency), and tOb = y^o 

confinement energy. Using (n, s, k) as the noninteracting eigenstate quantum numbers, where ft = (n, n') is the 
orbital Landau level index for the two decoupled ID simple harmonic oscillators, s = ±1/2 the electron spin eigenvalues 
along the total magnetic field, Bfot = (-B||,0,Bj^) and k the guiding center coordinate, one can obtain the following 
noninteracting eigenenergies, , for Hq^ of Eq. (1): 



n.s 



0Ji\n + - 



and the noninteracting orbital wavefunction: 



,/.y;,_,(r) = ^ ^^iiHcos^(x + llU) - sinez) ■ i^i',\smeix + l^k) + cos0z), 



(3) 



(4) 



where Ly is the system length in y direction and the function + l\k, z) has x and z components only; tan(2^) = 
—2uj\_uj\^l(uJ^ — w^). In Eq. (4), the function ijjn\x) is 



cxp 



'2lF 



(5) 



with li = ^J\lm*uji for i = 1,2, and Hn{x) is Hcrmitc polynomial. 

A typical noninteracting energy spectrum as a function of in-plane magnetic field, , is shown in Fig. 1 (with system 
parameters similar to the experimental values of Ref. [15]). In Fig, 1, two kinds of level crossing can be observed, one 
is in the weak B\\ (< 5 Tesla) region and the other in the strong B\\ (> 19 Tesla) region. In this paper we will denote 
the former to be an " intersubband level crossing" and the latter to be an " intrasubband level crossing" (see Table I) 
[43]. In this paper we adopt a convention that all single particle states in Eqs. (3) and (4) that have the same quantum 
number n (the first index of n) correspond to the same subband. With this definition of "subbands" , the characteristic 
energy separation of intersubband levels, Wi, is always larger than the energy separation of intrasubband levels, W2- 
It is useful to emphasize that the states of Eqs. (3)- (5) are the exact noninteracting eigenenergies for arhitrary values 
of in-plane field, perpendicular field, and confinement energy. In the second quantization notation, the noninteracting 
Hamiltonian can be written to be 



w 

— / . 

n,s,k 



w iv,t w 



n,s '"n,s,k'"n,s,k^ 



(6) 



where c^^^^(cS^s fe) creates (annihilates) an electron in the state (n, s, k). 

Before showing the full many-body interaction Hamiltonian, it is convenient to define the form functions [17], 



= I dx 



dze-"i^-e-'^'^^l{x - llqy/2, z)^'^^{x + llqy/2, z), 



(7) 



which are constructed from the noninteracting electron wavefunctions (Appendix B). In principle one could use 
the electron wavefunctions obtained either from a self-consistent Hartrec-Fock approximation [17,19] or from a self- 
consistent local density approximation [44,45] in Eq. (7) to calculate the form function. However, in Ref. [17] we 
have shown that the difference in the form function between using the noninteracting single electron wavefunctions 
and the self-consistent Hartrce-Fock wavefunctions is very small. Therefore we define the form function by using the 
noninteracting wavefunctions and consider only the diagonal part of the Hartree-Fock potential (i.e. the first order 
Hartree-Fock approximation) in calculating the ground state energy. In addition, we include the screening effect of 
positively charged donors by using the statically screened Coulomb interaction, V{q) = e^/eo(|g|^ -I- (27r/?s)^) for 
convenience, where is the effective screening length by donors outside the well and eo = 12.7 is the static lattice 
dielectric constant of GaAs. We believe this simple approximation should give reasonable results compared with the 
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more complicated calculation by including the donor density in a self-consistent approximation [44,45]. Using the 
notations defined above, the interaction Hamiltonian can be expressed as follows [17]: 

= — V V V (ff\p-iQ.{ki-k2)llJV,^ JV W,f W 

2f2 / ^ / ^ / y nin2,n3n4,\'i/ ^ ^ni,<Ti,ki+qy/2^n2,CTi,ki-qy/2'^n3,(T2,k2-qy/2^ni,a2,k2+qy/2 
ni,---,n4 CTi,a-2 ki,k2,q 

= ^ V V V Vi^^ ^ ^ (ff,]p-i'lAki-k2)llJV,i „W W,\ W /ci\ 

2Qj^ / ^ / ^ / J ^ nin2,n3n4,Vi-^ I ^ ni,(Ti,ki+qy/2^n2,cri,ki-qy/2'^n3,a2,k2-qy/2^n4,c72,k2+qy/2J X'-'J 

ni,---,n4 cri,(T2 /ci,fc2,Q_L 

Where Vif^^ .3^,(9) = ^(«X„^,(-gX„^,(a) and V;]f^„^3^,(gl) ^ i"^ E,. V;rn2,n3n4(«)- " = ^^^y^^ the 
usual normalization volume and Ox = L^Ly is the 2D normalization area. If we choose the alternate Landau gauge, 
A[2,](r ), the expression of will have the phase factor, e~"'»('^i~'^2)'oj rather than e^'^^'^^^^^^'^^o above. 

For the convenience of later discussion, we can express the Hamiltonians shown in Eq. (6) and Eq. (8) as follows: 

-f^t^t = X^^S,r^^,nCT(0) + 2J2~X! ^ X] ^*5rn2,S3"4(^L)Ps^<7i,n2CTi(-^L)p^a2,n4<T2(9-L)' (9) 
n,(T q±_ ni,-",n4 fiiCr2 

where 

izf \ = Sr^ Jqu,kllJV,\ W /ir>\ 

yniCr^,njCr2Vl-^) — ^ni,ai,k-qy ,a2,k+qy /2 

k 

is the density operator for the wide well system. 



B. Double Well System 

For a double well (bilayer) system, we assume both wells are of zero width in their growth direction so that the 

in-plane magnetic field does not change the electron orbital wavcfunctions as in wide well systems. For most bilayer 
problems of physical interest, neglecting the indivitual layer width is an extremely reasonable approximation. In the 
presence of an in-plane (parallel) magnetic field, the; tmnieling amplitude in the Hamiltonian acquires an Aharonov- 
Bohm phase factor to satisfy the gauge invariance of the whole system [22,36,27]. This, in fact, is the main efi^ect of 
the in-plane field in the double well system for our purpose. To describe this important effect, we first express the 
Hamiltonians of a double well system in a special gauge [26]: A'^y^{f) = {0,B±x,Bxy — Byx), where and By are 

the X and y components of the in-plane field (i.e. _B|| = (Ex, By,0)), by using a conventional basis of electron states, 
(n, /, s, k), where n is Landau level index in a single layer, I = ± is the layer index for right(left) layer (sometimes also 
called up/down layers in analogy with electron spin), and s and k are the same as before. The Aharonov-Bohm phase 

factor, exp[— z4>g I±d/2 ^[y] z(^J-' ^)'^^] (where d is layer separation and <I>o — hc/e is fundamental flux quantum), is 
then introduced for electron tunneling from one layer to the other. In the bilayer Hamiltonian we need to keep the 
tunneling term only for the electrons in the highest flUed Landau level, and therefore the noninteracting Hamiltonian 
in the second quantization representation becomes: 





^\ 1 











D,] D 
^n,l,s,k^n^l,s^k 



nj,k,s 

-rjV,P2^|e '^Ar, + l,s,fc-Px/2''JV,-l,s,fe-|-P„/2 "I" ^ '-N-l,s,k+P^/2^N,+l,s,k-P^/2 ' \^^) 



s.k 



where c^'/g ^(c^; « k) ^ the electron creation (annihilation) operator of state (n, I, s, k), described by Ly ^^^e''''''^ipn\x+ 
klQ)S{z =F d/2), where 'ipn\x) is the same as Eq. (5) with k replaced by magnetic length ^o- -Pi = {Px,Py) = 
27rdi?||/<I'o is a characteristic wavevector introduced by the in-plane field. The effective tunneling amplitude t^.p is 
f e~^^'o/"'L^(P^^o/2) with t being a parameter for interlayer tunneling, where P = \P±\ and L^{x) is the generalized 
Laguerre polynomial. Following the convention in most of the existing literature, we will treat tjv,p as a whole to 
be an independent variable in this paper and hence neglect all its level and magnetic field dependence for simplicity 
[46,47]. 

Similar to Eq. (7) we define the following form function for the double well system: 
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where (l)Zkir±) ^ Ly^^^e'''y^/n \x + kll) is the one electron wavefunction of Landau level n in a single well. The 



z 



component of the wavefunction can be integrated first and the resulting explicit formula of y4.^.^ .(gj_) is shown in 
Appendix B. The interaction Hamiltonian in this basis is then 

rD = ^ V V (ff,\p-iq.{ki-k^)llD,\ D D,\ D 

2r2j_ ^ Z-^ ''nin2,ra3"4V^-'-/ ni,ii,si,fei+g^/2^n2,h,si,fei-9a/2^ri3,;2,S2,fe2-ga/2 "4,i2,S2,fe2+g!,/2 

m-'-m li,l2 si,S2 ki,k2,q± 

" ^J^X]X]^n^n2!"3"4(iL)/>?i;i,„2ii(-^L)P^3'i2,„4i2(^L) (12) 
q± hM 

where y„^;^i;5,23„^(gj_) = V{qi_)A^^^^{-q^)A^^^^{q^)ex.^{-\\h-l2\d\qi_\), including the contribution of the z compo- 
nent of the electron wavefunction; p^^i^^iA^^) = T,s,k ^''^''''^°%.uhs,k-qj2'^2,i,s,k+qy/2 is ^'^^ electron density operator 
of layer I and spin s. 

Sometimes it is more convenient to rewrite the Hamiltonian in terms of the eigenstates of the noninteracting 
Hamiltonian, which are defined as: 

^ gi(fe+p,/2)p„;^/2 ^ 

^n,+l,s,fe — \^n,+l,s,k+P^/2 + °'n,-l,s,fe+P«/2 j 

'^n,-l,s,k ~ \^n, + l,s,k-P^/2 ~ ^n,-l,s,k-P^/2 J ' (1^) 

where a^^g ^ is the new electrons creation operator of a symmetric (a = +1) or an antisymmetric (a = —1) nonin- 
teracting state, because they are also the eigenstates of parity (space inversion) tranformation (cj^ ^ ^ 



k ^ ^n. — l,s. — k J ' 

Note that the description in terms of the left/right layer index / (in Eq. (11)) or the symmetric/ antisymmetric state 
a are completely equivalent since they are simple linear combinations of each other (cf. Eq. (14)). 
In the a basis, the noninteracting Hamiltonian of Eq. (11) is 

= X/ ^n,a,sC'li,a,s,k^n,a,s,k^ (15) 
n,a,k,s 

where the noninteracting electron eigenenergy, -B^'^g, is 

En:Zs - (" + 1/2)^^^ - aAsAs/2 - su,, (16) 
and AsAS = 2t is the tunneling energy. The interacting Hamiltonian in Eq. (13) then becomes 



al,a2,q± m ■•■"4 



A=±l 



(17) 



where the new density operator, p,^J°|(9l) = T,s,k e''"''''°«ii,c«i,s,fc-,,/2«n2,a2,s,fe+g^/2 and the three interaction matrix 
elements in Eq. (17) are as follows: 

V^fn2,nsnAQ^) ^ \ (v;'?;t!n3n4 (^l) ± ^Os{{q^ ■ P±)ll)K'',n\:n,nM)) (18) 
Kn2,n,nM^) ^ ^ sin((gl • PMV^^t^n.nA'l^) ■ (19) 

As will become clear from the later discussion in this paper the basis described by Eqs. (15)-(19)) is often more 
convenient for constructing trial wavefunctions, especially for an even filling system in the presence of in-plane magnetic 
field. 
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C. Symmetry properties of the systems 



We now discuss the symmetry properties of the wide well systems and the double well systems in the presence of 
in-plane field. Wc first consider parity (full space inversion) symmetry, translational symmetry and spin rotational 
symmetry individually, which exist in both the systems of interest (i.e. wide well and double well). Then we discuss 
the ^;-parity (reflection about x — y plane) and the isospin rotation symmetries, which exist only in the double well 
system in certain situations. 

In this paper, we define the parity operator, V, so that it reverses the full space coordinates from r = {x,y,z) to 
— r = {—X, —y, —z). Therefore, in the second quantization notation, we have 

-PcZ^k-p-' = (20) 
for the wide well system (see Eq. (4)), and 

'P'^n,l,s,k'P~' = °' 'Pan,a,s,k'P-' = (-l)"aa„,„,«,_fc (21) 

for the double well system in the two different basis. It is easy to sec that the Hamiltonians of these two systems 
shown above are not changed under the parity transformation by using the identities (sec Eq. (Bl)): AW^-,^^{-q) = 
(-l)"i+"'i (-l)"2+r4^M^_^^ (- j^D^^^ i-q±) = (-l)"'+"'^^'i„2 (^l)- This simple result is vahd even in the presence 
of the tilted magnetic field, if only we consider a symmetrically confined (but not necessary parabolic) potential well in 
the wide well system and no external bias voltage in the double well systems. Thus, parity is rather general symmetry 
property of the physical systems under consideration. 

To consider the 2D translational symmetry in the x — y plane, it is more convenient to consider the many-body 
system in the first quantization representation, and introduce the total "momentum" operator, M. [48,17]: 

M = ^[p^ + ^A{n) - ^Btot X fi] , (22) 

i 

where pi and fi are the momentum and position operators of the ith electron. A is the vector potential and Btg^ = 

(i3|| , 0, B±) is the total magnetic field. Then it is straightforward to see that the x and y components of A4 commute 

with the Hamiltonian, H = J2iiPi + e^(^i)/c)^ + J2i Uc{zi) + \ Yliijtj ^{^i ~ ^j)- This result is true even in the 
presence of an in-plane magnetic field and for arbitrary shape of the electron confinement potential in z direction. 
Therefore it applies equally well to wide well systems and double well systems. We can then define the 2D translation 
operators in the x and y directions as 

%{R^) = exp {-iR^M^ , (23) 
fy{Ry) = exp (-iRyMy) , (24) 

which form closed translational symmetry group individually [49] (they do not commute with each other due to the 
presence of the magnetic field). We first study how %;^y{Rx,y) can transform a noninteracting electron eigenstate of a 
wide well system, (f^i;{r) = ^^^e''^^$^(a; + l^k, z), with a shift, = {R^, Ry), in the x — y plane. In the Landau 
gauge, A[yj{f)), we obtain 

TM^'SAn = e-^^«^/'°.^^,fe(a; - R,,y,z) = cPZ,k-R./ilin, (25) 
fy{Ry)^JAf) = ^'Y.ix, y - Ry, z) = e-'^^ycf^,k{r). (26) 
In the second quantization representation, it is equivalent to 

fy{Ry)cl^,^A{Ryr' = e-^'^^^'cl,,,,. (28) 

In other words, T^XR,.) shifts the guiding center coordinate and Ty(Ry) adds an additional phase factor of the 
electron operators respectively. As discussed in Section II B, in the case of double well systems we use the gauge 
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j4|^j(f) = {0,B±x,Bxy — Byx) (see the paragraph above Eq. (11)). Therefore the "momentum" operator in this 
gauge becomes M. = {Px + v/Iq ~ ^Py/d^Py + zPx/d,Pz) and the translational operators are 



%{Rx)<l>^,,{f^)S{z - ld/2) = e^'^«^''/V^,,_«^/;.(r-l)<5(^ - ld/2), (29) 

t,(i?,Xfe(r-l)<5(^ - ldl2) = e-''«^^»/^e-''=«''<,(r-l)(5(^ - ld/2). (30) 
In the second quantization representation, they are equivalent to 

t(i?.)c^,,..fct(i?.)-^ = e^'^^^''/^<,,,,+^^/,., (31) 

%{RyKi,sAiRyr' = e-^'^^-'^^/'e-'^^^^,,,,. (32) 



Comparing Eqs. (31)-(32) to Eqs. (27)-(28), we find that an additional phase factor appears in the layer index basis, 
which is related to the Aharonov-Bohm phase factor and hence confirms that the tunneling Hamiltnonain in Eq. (11) 
commutes with the translational operator, If we use Eq. (14) to change the above results (Eqs. (31)-(32)) to 

the symmetric-antisymmetric basis, we obtain 

Tx{Rx)an,a,s,k%:{Rx) ^ = fln^Q.s.fe+iJ^//^ , (33) 
y^Ry^^^^ocjSyk'^yi^Ry) =6 ^^^n,a,s,fc? (^^) 

which are the same as Eqs. (27)- (28) obtained by replacing ^ j, '^y cin a s k- 

As for the spin rotational symmetry about the total magetic field direction with an angle x, we can simply apply 
the following operator on the spin wavefunction, 

U{x)=e-'^'% (35) 

where Sz is the z component of spin operator with the spin z axis being along the direction of total magnetic field. 
Since our noninteracting eigenstate is always the eigenstate of Sz, the spin rotational operator, li{x), just gives 
an additional prcfactor, e^^^''^ for spin up(down) electron annihilation operators. The spin rotational symmetry is 
conserved by the Hamiltonians that we consider. 

Now we discuss the symmetry properties that exist only in the double well system but not in the wide well 
system. First, if there is no in-plane magnetic field (i?|| = 0), the double well system has two separate parity 

symmetries: one is the z-parity symmetry, Vz (changing z to —z) and the other one is the in-plane parity symmetry, 
Pxy (changing {x,y) to {—x,—y)). In the second quantization representation, we have PzC^i s k'^z^ = -i s k ^^'^ 

VxyCni s k'^xy = (~l)"c^i s -fc- is casy to show that both Vz and Vxy commute with the total Hamiltonian of the 
double well system in the absence of in-plane magnetic field. When an in-plane magnetic field is applied, however, 
none of them is conserved due to the Aharonov-Bohm phase shift in the tunneling amplitude (see Eq. (11)), while 
their product, the full space inversion V = Vz'Pxy, is still conserved as discussed above. Second, if there is no electron 
interlayer tunneling (A^^g = 0) in the double well system, the electron number in each layer is also a constant 
of motion. In the isospin language of the layer index basis, such conservation is described by an isospin rotational 
symmetry about the isospin z axis, UisoiO = ^^P[~'''i{(^n'+,s,k^n,+,s,k~^n,'-,s,k'^n,-,s,k)/'^]- The spontaneous interlayer 
coherence will break such continuous symmetry and give rise to a Goldstone mode. Such U{1) symmetry breaking 
has been extensively studied in the literature (see Ref. [22] and references therein), and hence wo will not seriously 
address this issue in this paper. However, since our interest in this paper is to study the charge and spin density 
wave phases induced by the in-plane magnetic field, we will always consider systems with finite tunneling and finite 
in-plane magnetic field throughout, so that the above unique symmetries (Pz, Pxy, and Uiso) actually do not exist 
in the double well systems we consider in this paper. As a result, both the wide well systems and the double well 
systems we study will have full space parity symmetry (V), two-dimensional translational symmetry (Tx^y), and spin 
rotational symmetry about the total magnetic field (U). We will then discuss how a stablized many-body coherent 
wavefunction can break these symmetries near the level crossing (or level near degeneracy) regions. 

III. VARIATIONAL WAVE FUNCTIONS 

In this section, we propose trial many-body wavefunctions to variationally minimize the Hartree-Fock energy of the 
integer quantum Hall systems in the level degeneracy region: a wide parabolic well and a double well system at both 
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even and odd filling factors. Stripe phases, that we are primarily interested in, have translational symmetry along 
the longitudinal direction of the stripes. They are conveniently discussed in the Landau gauge, for which particle 
momentum is a good quantum number in the direction of the stripe. In this paper we only consider stripes that 
are parallel or perpendicular to the direction of the in-plane magnetic field. Therefore for the in-plane field in the 
X direction we choose gauges A[yj{f ) = {0,B±x — B^\z,0) and A^x]{f ) = {~B±y,—B^\z,0) to describe phases with 
stripes along y and x axes respectively. By comparing the HF energies of these two stripes phases (one is along y 
direction obtained by using gauge A[y^{f), and the other one is along x direction obtained by using gauge A[j.](r)) 
and the HF energies of other nonstripe phases (discussed later), we can determine if a stripe phase can be stabilized 
and which direction is energetically more favorable for stripe formation. If no stripe phase is energetically favorable, 
the nonstripe phases obtained in these two gauges are identical (as they should be due to symmetry). In the rest of 
this paper, for the sake of brevity, we will show equations and formulae for the trial wavefunction and the related HF 
energy of a wide well system only for stripes along y axis (perpendicular to the in-plane field) in the gauge A^y]{f), 
although (we emphasize that) we always consider both stripe directions in our variational calculations. 

For a zero width double well (bilayer) system, the trial wavefunction for a stripe phase along the in-plane field 
direction {x) can be studied in a much easier way. Using the fact that the in-plane magnetic field does not change the 
electron orbital wavefunctions in the zero width wells, we can simply rotate the in-plane magnetic field, By, from the 
direction along x axis to the direction along y, keeping all the trial wavefunction obtained by the conventional gauge, 
j4[y](r ), the same. We then equivalently obtain the HF energies of a stripe phase with a conserved momentum either 
perpendicular or parallel to the in-plane magnetic field. Their energy difference results from the Aharonov-Bohm 
phase factor in the tunneling amplitude, and that is the reason we keep both and Py in the bilayer Hamiltonian 
in Eq. (11). 

In the rest of this section, we first show the trial wavefunction for the two level crossing (or level near degeneracy) 
situations in Section III A. We also consider several special many-body phases it generates and study their physical 
properties. In Section IIIB we propose wavefunctions for the four level degeneracy, involving the second highest filled 
level and the second lowest empty level about the Fermi energy for the double well system aX v — AN + 2 and for 
the wide well system &t v ~ 2N + 2, leaving out all other lower filled levels and higher empty levels as irrelevant core 
states. In Section HI C we develop a perturbation technique to investigate the possible instability toward the stripe 
formation from a uniform coherent phase. We can use this method to study the existence of a stripe phase and its 
basic properties without optimizing the whole HF energy if the stripe formation is a second order phase transition. 



A. Two level degeneracy: variational wavefunctions for odd filling factors 



1. Wavefunction 



We propose the following trial many-body wavefunction to incorporate both isospin stripe and isospin spiral orders 
simultaneously: 



e''=««'o/2-Fi7/2cos(^fc/2) e- 
-e'*^'3-'o/2+i7/2sin(V;fe/2) e" 



ikQJ„/2-ij/2 sin(^fe/2) 

ikQa,il/2-i'r/2 cos(V;fe/2) 



'-it,k-Qy/2 
_ %,k+Qy/2 



(36) 



where we use f (JJ.) to denote the isospin up(down) state, which can represent subband level, spin and/or layer indices 
depending on the systems we are considering (see Table I). The spiral winding wavevector, Q±, stripe phase function, 
ipk, and the additional phase, 7, are the variational parameters to minimize the total Hartree-Fock energy. If the 
system has isospin rotation symmetry around the isospin z axis (e.g. double well systems a,t u = AN+l in the absence 
of tunneling), the new ground state energy obtained by the trial wavefunction above will be independent of 7. If the 
system has no such isospin rotation symmetry (e.g. double well systems at 1/ = AN -|- 1 in the presence of tunneling 
or the wide well systems a.t v = 2N + 1), the new ground state energy will then depend on 7, selecting some specific 
values of 7 to minimize the total variational energy. We will discuss both of there cases in details in our HF analysis in 
Sections IV- VII. In general this two level coherence approximation should be good for a system near the degeneracy 
point, especially for an odd filling factor. For systems with even filling factor, the inclusion of the next nearest two 
levels about Fermi energy becomes crucial as will be discussed later. 
For the convenience of later discussion, we define following functions: 
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^^"^ fe 
't' k 

Qsiq.) = ^Y1 e'"'"'^ sin(V'fe/2) cos(^fe/2), (37) 

Note that when tf^k = ■^o as a constant, 0i(gx) = ©i(0)(5g^,o (for i = 1,2,3) and ©3(0)^ = ©i(0)O2(0). The physical 
meaning of these functions are the following: 0i(O) describes the density of isospin-up electrons, ©2(0) describes the 
density of isospin-down electrons, and ©3(0) measures the coherence between isospin up and isospin down electrons. 
For the periodic function, ipk, without any loss of generality, we can assume that it is a real and even function of the 
guiding center coordinate k so that Qi{qx) are all real quantities. 

2. Isospin phases and their physical properties 

To understand the physical properties of the trial wavefunction proposed in Eq. (36), we first define the following 
generalized density operators (isospin index, J = ±1 =-fl- (JJ.)): 

Pl.lM) = E^'"'^'°4„fe-,,/2C7.,fe+,,/2, (38) 
fe 

which generalizes Eqs. (10), (13) and (17)). The isospin operator at the in-plane momentum, q±, can be defined as 
i{q±)^l J dz J df^e-''---- E 4.,kAui.^i.MK,kAn*KMi^) 

h.h fei,fe2 

= E '^^1,^2^/1/2 (91)^/1/2(51,0), (39) 

I1J2 

where aj^ are the Paidi matrix elements, and Aj^j^{q) is the generalized form function as has been specifically 
defined in Eqs. (7) and (12) for wide well systems and double well systems respectively. 

Using the trial wavefunction, |^'g)j in Eq. (36), we obtain the following expectation value of isospin components: 

(4(5*)) = I [A^f{q±,0){^G\pni'l±)\'^G) - A4^((f±,0)(*G|PM(9l)l*G)] 

At^(5-l,0)e-*«==Q'''o/2ei(fe) - ^M(3l,0)e'«^«'''c/2©2(fe)l (40) 



and 



(X+(ql)) - Ai^4(gl,0)(vl'G|m('?l)l*G) = A^A„.Q.^M(9l,0)©3(g. - Q.)e-"^. (41) 



There are two classes of interesting phases one can find from Eqs. (40) and (41). First, if tpk = V'o is a constant, 
the z component of the mean value of isospin, {Iz{q±)) oc (^9^,0^9^,0, is uniform over the 2D well plane, while its 
transverse components, {T±{q±)) oc Sq^^g^Sq^^Q^, select a certain wavevector, Q± = {Qx,Qy), in the guiding center 
coordinate for winding, i.e. the isospin x and y components oscillate in the real space at wavevector, Q± (see Fig. 
3(a)). Therefore we define Q± as the wavevector of isospin spiral order. Secondly if ipk is a periodic function of 
k, the isospin polarization, (^^(r _L)), will oscillate along x direction (but uniform in y direction) according to Eq. 
(40), characterizing a stripe phase with normal vector n\\x. Choosing the other gauge, A[^^(f), in which the particle 
momentum is conserved along x direction, we can construct a stripe phase along x direction and isospin (^^(r _L)) 
then modulates in y direction as shown in Appendix A. 

It is also interesting to investigate the local charge density distribution in the top Landau level, P;ocai(^J-)> ^y using 
the trial wavefunction, I'^g), in Eq. (36). It has contributions from both the isospin up and isospin down electrons: 
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Plocalir±) = -^Yl e'^"^'^"^ {^G\A^fiq±, 0)pui^±) + ^u(9±,0)Pu(^±)\^g) 

= - ^sfe) [^u^Qx, 0, 0) cos{q^{x - Qyll/2)) - A^^{q^, 0, 0) cos{(i^{x + Qyll/2))\ , 

= po + Pex{r±) (42) 

where N^/^±_ = {2'kIq)~^ is the average electron density in each Landau leveL Eq. (42) shows that if Q2{Qx) selects 
a specific wavevector (i.e. Vfe oscillates periodically with a characteristic wavevector) at Qx = Qn = ^nn/a, where a is 
the period of stripe, the extra charge density, Pex{r±), will be nonzero and a periodic function in the real space. In 
Appendix C we will show that this extra change density is related to the charge density induced by the topological 
isospin density by generalizing the theory of skyrmion excitations developed in Ref. [18] for a double well system. 
Therefore in the rest of this paper we will denote such a CDW state as a "skyrmion stripe" phase. In general we 
should have skyrmion stripe phases whenever the stripe normal vector n (i.e. the modulation direction of (Xz)) is 
perpendicular to the spiral wavevector Q± [8]. 

According to the above analysis, we can consider the following six different phases obtained from Eq. (36): when 
ipk is a constant, the wavefunction of Eq. (36) can describe three nonstripe phases: (i) fully (un)polarized uniform 
quantum Hall phase for tjjk = (0)7r (with Q± being arbitrary), (ii) coherent phase for Q± = and ipk 0,7r, and (iii) 
spiral phase for finite Q± and ipk 7^ 0, n. When tpk changes periodically with k, three kinds of stripe phases arises: (i) 
when Q± is perpendicular to the stripe normal direction, we have a skyrmion stripe described above; (ii) when Q± is 
parallel to the stripe normal direction, we have a spiral stripe phase, and (iii) when Q± = we have a coherent stripe 
phase. For simplicity, in this paper we will only consider stripe and spiral orders (if they exist) with characteristic 
wavevector n and Q± being either perpendicular or parallel to the in-plane magnetic field, which is fixed to be in 
the positive x axis. A simple visualized cartoon of the spiral winding and the three stripe phases is shown in Fig. 4. 
We calculate the HF energy given by the six different phases generated from the trial wavefunction, Eq. (36), and 
compare them to get the minimum energy for the true ground state. 



3. Symmetry properties of the trial wavefunction 

In Section II C we have shown that the systems we are considering have three kinds of symmetries: parity (full space 
inversion) symmetry, two-dimensional translational symmetry, and spin rotational symmetry. Here we show how a 
general isospin trial wavefunction in Eq. (36) is transformed imder these symmetry operators defined in Section II C. 
This will enable us to identify the broken symmetries in each of the isospin phases discussed above. 

First it is easy to see that the conventional integer QH state (i.e. isospin polarized with tpk = or tt), |^o) = 
rifc '^1(4) fcl-^-^)' eigenstate of all the three symmetry operators {P, %:(y), and U), by applying the equations in 
Section II C directly. (For simplicity, in our discussion below, the isospin basis of double well system a,t v = AN + 1 is 
chosen to be the symmetric-antisymmetric basis, i.e. the noninteracting eigenstate basis. The symmetry properties 
of the isospin coherent phases constructed in the layer index basis will be discussed later in Section VI D). Therefore, 
the isospin polarized state, I^I/o), does not break any symmetry properties as expected. 

Now we consider a general many-body state described by the wavefunction Eq. (36) with variational parameters 
obtained from minimizing the HF energy: ip^ = tpl 0,7r; Q± = Q*j_, and 7 = 7*. Applying the parity operator, V, 
translation operator, Tx(y), and spin rotation operator, U, on |^'g(V'o, Ol, 7*))) we can obtain respectively 

Vl^cm, Ql,!*)) = (-i)^1*G(V-fe, -Ql, 7* + tt)), (43) 

t„(ii,)|*G(V'fc,Ql,7*)) = \^G{rk_f,.n^/i2,Ql,7*-QaRa)), (44) 

ZY(x)|*G(V'fc, Ql, 7*)) = l^cCV-fc, Ql, 7* - X{si - S2))), (45) 

where a = x or y, R± = {Rx,Ry), Si(2) = ±1/2 is the spin quantum number of isospin up(down) state, and n is 
the stripe oscillation direction (it is perpendicular to the direction of stripes, e.g. fi = ^ for stripes along y direction 
described by the Landau gauge, A[y]{f)). In Eq. (43) we have assumed that the signs for the isospin up state and 
isospin down state are opposite after parity operation (we will show that this is always true in the level crossings 
considered in this paper). According to Eqs. (43)- (45), we find that only parity symmetry is broken in a coherent 
phase (V'fc = V'o 7^ 0> ^ti" and Q*^ =0), if Si = S2- Spin rotation symmetry is also broken if the spin quantum number 
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of the two crossing levels (isospin up and down) are different (si ^ S2), which is true only for the level crossings in the 
even filling systems. When we consider the spiral phase with ip'^ = ip^ 0,tt and Q*j_ ^ 0, we find that in addition 
to the broken symmetries discussed above (parity symmetry and spin rotational symmetry if si ^ S2), it breaks 
translational symmetry in the direction of isospin winding (i.e. along Ql)- From Eq. (44) we immediately see that 
the wavefunction has a period 2tt/\Q*^_\. For the coherent stripe phase {tjjl modulated periodically and Q*^ = 0) we 
have broken parity, spin rotational symmetry if si ^ S2, and translational symmetry in the direction of n. For a spiral 
stripe, the spiral winding direction is parallel to the stripe oscillation direction (QlUn), and therefore the translational 
symmetry is broken only in one direction, while a skyrmion stripe (Q*j_ _L n) breaks translational symmetries in both 
X and y directions. Both the spiral stripe and the skyrmion stripe also break the parity symmetry and spin rotational 
symmetry if si ^ S2. 

A fundamental quantum mechanical principle stipulates that when a system undergoes a quantum phase transition 
to break a symmetry (i.e. the ground state wavefunction is not an oigcnstate of the symmetry operator), the new 
symmetry-broken ground state will have additional degeneracy associated with the spontaneously broken symmetry. 
This result is also obtained in our HF energy calculation shown later. In Table I we list the many-body states obtained 
by our HF variational calculation and the resulting broken symmetries. We will discuss each of them individually in 
the following sections for different systems and then compare the results with each other in Section VIII. 

B. Four levels near a degeneracy: variational wavefunctions for even filling factors 

1. Wavefunction for double well systems 

For a double well system at v = AN + 2, the complete filled core levels are the lowest 4A^ levels, while the top two 
filled levels may coherently hybridize in some situations with the empty levels above the Fermi energy (we consider 
only the two lowest empty Landau levels in the context motivated by the scenario originally discussed in Ref. [33]). 
We label the states involved in forming a many-body state as follows: the second highest filled level is denoted to 
be level 3, the highest (top) filled level is level 1, the lowest empty level above the Fermi energy is level 2, and the 
second lowest empty level is level 4 (see Fig. 5(a)), and the other higher empty levels are assumed to be irrelevant. 
Naively one may think that we may construct a trial wavefunction similar to Eq. (36) to consider interlevel coherence 
only between level 1 (the highest filled level) and level 2 (the lowest empty level), because the single electron energy 
separation between these two levels is the smallest energy scale near level crossing. However, it was explicitly shown 
that in a double well system a.t u = AN -\- 2, the contribution of the second highest filled level (level 3) and the second 
lowest empty level (level 4) could make crucial contribution to the coherent hybridization between level 1 and level 
2, resulting a novel canted antiferromagnetic phase (CAE) with broken spin symmetry [33]. This symmetry-broken 
coherent phase cannot be obtained if one considers the coherence between the most degenerate pair only (i.e. level 1 
and level 2) via Eq. (36). Another way to understand why it is natural to consider four rather than two levels around 
the Fermi energy in creating a many-body state for the double well system aX u = AN -I- 2 is to note that these four 
levels are separated from the other ones by a large cyclotron energy, whereas they are separated from each other only 
by (much) smaller energies of Zeeman and tunneling splittings. Therefore in this section we will consider the mixing 
of the four (rather than two) levels closest to the Fermi energy to study the trial wavefunctions for the double well 
y = AN + 2 system near the level crossing region in the presence of the in-plane field. 

Let us first consider the uniform ground state of a double well system at v = AN + 2 without any in-plane magnetic 
field. We assume that the cyclotron resonance energy, is much larger than the tunneling energy, Ag^s, and the 
Zeeman (spin-splitting) energy, lOz , so that the two highest filled levels and the two lowest empty levels have the same 
orbital Landau level and all other levels can be treated as incoherent core states not actually participating in the 
level-crossing hybridization process. Single particle states for noninteracting electrons are shown in Fig. 5(a) in the 
case when the tunneling energy, Ag^*?, is larger than the Zeeman energy, oj^- We use (a, s) to label the four levels 
around Fermi level under consideration, where a = ±1 is the quantum number of parity symmetry (which in this case 
is the reflection symmetry about x — y plane, z —z, or the interchange of layers, see Section II C), and .s =| (j) is 
the spin quantum number. For convenience of comparison and later discussion, we define level 1 ((a, s) = (+, |)) and 
level 2 ((a, s) = (— , \)) to be the isospin up and isospin down state respectively as shown in Table I, because these two 
levels, being closest to the Fermi energy, are obviously the most energetically relevant ones compared to the other two 
levels in this four level scenario (see Fig. 2(f)). Many-body states that appear in this system correspond to mixing 
some of these single particle states and describe the breaking of certain symmetries in this problem. As discussed by 
Das Sarma et. al. [33] in the context of the canted antiferromagnetic state in bilayer systems, the states that are most 
likely to be hybridized by Coulomb interaction are the ones that are closest to each other for the noninteracting system: 
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state 1 = (+, J,), and 2 ~ (— , "f). When the expectation value (cj ^c^ ^ is finite and independent of k, the system has 
canted antiferromagnetic spin order. In the two layers the transverse components of the spin point in the opposite 
directions. Such an order parameter breaks the spin rotational symmetry around z-axis and and the parity symmetry 
V = VzVxy (We use the full parity symmetry, V, rather than Vz, since the latter is not conserved when we include 
an in-plane magnetic field later, see Section II C.) More precisely, under the spin rotation, U{x) (defined in Eq. (35)), 
we have U{x)a\ }fi~^{x) = e~'^a| f^, and under parity transformation, we have "Poj kV~^ = —a\ _^a2 _f.- 
So the order parameter has spin 5^ = 1 and is odd under parity. However, the operator 04 f.a^ j. has exactly the 

same symmetry properties a\ ^a^ so in a canted antiferromagnetic phase both of them acquire finite expectation 
values. Therefore it would be insufficient to consider mixing of the states 1 and 2 only and treat levels 3 and 4 as 
frozen [33] in the double well system at even filling factors. On the other hand hybridization between any other pair 
of levels does not take place since the appropriate expectation values would have symmetry properties different from 
the transverse CAF Neel order. For example a\a^ has the right parity symmetry but wrong spin symmetry, and a\a^ 
has the correct spin symmetry but wrong parity. As was discussed earlier [33-35] and as we will demonstrate below, 
the physical origin of the CAF phase is the lowering of the exchange energy due to the additional spin correlations 
associated with a\ f.a^ ^, present above. 

To construct an appropriate trial wavefunction in the presence 0/ in-plane field, we first consider a uniform phase. 
If we make the assumption that the order parameter breaks the same symmetries as in the case = 0, it has to 
have spin 5^ = 1 and be odd under parity. The operators, a^__,^ ^. defined in Eq. (14), have the same transformation 

properties under U{x) and V (see Section II C) with or without B|| . The identical nature of broken symmetries in the 
cases B|| = and B|| ^ implies that a\ ^.Oj and a\ /.a^ ^ acquire expectation values even in the presence of a finite 
in-plane magnetic field. Therefore we can still apply the ansatz of Ref. [33] and consider the hybridization of level 1 
with level 2 and level 3 with level 4 separately to construct a trial wavefunction wavefunction for the uniform phase 
in the presence of in-plane field. In addition, as suggested by the wavefunction of odd filling system in Eq. (36), this 
hybridization may be further extended to a slightly more complicated form to reflect the possible uniform winding of 
the transverse Neel order [36] and/or the stripe order to break the translational symmetry. We then propose a trial 
wavefunction for a double well system at z/ = AN -|- 2 in the presence of in-plane magnetic field as follows: 



,fe"2, + ,fc 



\LL) 



2,-|-,fe 

r.t 



^^kQJl/2+^^/2 cos(7/;fc/2) e-*fcQ.io/2-n/2 sin(V;fe/2) 

_pikQ.ll/2+if/2 sin(V)fe/2) e-*'=Q-'o/2-i7/2 cos(Vfc/2) 







^N,+l,l,k-Qy/2 
k+Qy/2 
^N,+l,t,k-Q'/2 
^N-l,l,k+Ql,/2 






^ikQ'jl/2+iY /2 cos{i}'j2) e-*'3»'o/2-i772 sm{il)'j2) 

_^ikQ'jl/2+W /2 sm{^j;'j2) e-'*='3x'o/2-i772 cos(V'fc/2) 



(46) 



where -0^, Qj_ = {Q'xiQy)i and 7' are four additional parameters to be determined variationally. We allow V'fc and ip'j^ 
to be arbitrary periodic functions of fc, in order to consider the possible stripe formation. Note that the original 4x4 
matrix representation of an unitary transformation has been reduced to an effective block-diagonalized matrix form, 
and reduces to the uniform wavefunction originally proposed in Ref. [33] if we put Q± = and ijjk = ^Oi a constant. 
The fact that the uniform wavefunction {Q± = 0, ipi. = tpo) turns out to be an excellent description [33-35] for the 
corresponding B|| = case leads us to believe that the variational symmetry-broken wavefunction defined by Eq. (46) 
should be a reasonable generalization to study the many-body phases in the presence of an in-plane magnetic field, 
for the = AN + 2 bilayer system. 

The specific symmetry-broken form of the wavefunction in Eq. (46) allows us to introduce the concept of " double 
isospinors" to describe the coherence in the four levels near degeneracy, because each degenerate pair (i.e. levels 1 
and 2, and levels 3 and 4) form two distinct isospinors in the state defined by Eq. (46). The exchange energy between 
the two isospinors may stabilize a symmetry-broken phase, and it reaches its maximum value if the two isospinors 
have the same stripe period but have opposite spiral winding wavevectors. (This result is explicitly obtained later in 
our numerical calculations of Section VII.) The symmetry properties of this four level mixing coherent wavefunction 
(Eq. (46)) are identical to those discussed in Section III A 3, and we do not discuss it further here. 
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To develop a deeper understanding of the trial wavefunction proposed in Eq. (46), it is instructive to transform 
Eq. (46) back into the layer index basis, which, while not being a noninteracting energy eigenbasis, is physically more 
appealing and easier to visualize conceptually. For the convenience of comparison, we let Px = Qy = Qy = and 
7 = 7' = 0, i.e. the in-plane magnetic field chosen to be in the y direction. (Here we have used a known result that 
the winding vector is always perpendicular to the in-plane magnetic field in the double well system, which is justified 
by the numerical results shown later in this paper.) Combining Eq. (14) and Eq. (46) we obtain the ground state 
wavefunction as follows (the Landau level index, A'' is omitted): 

= n n \LLh (47) 

i=l,2 k \l,s J 

where z£l ^ = cos(V'fe/2), ^ = — sin(i/'fe/2), z'f^^ ^ = ± sin(^/'^/2), and z'f^^ ^ = cos(i/'fe/2), and their phases are 
respectively 

{TPy + Qx)/2 
{TPy - Qx)/2 

(TPy - 

iTPy + Q',)/2. (48) 

One can see that the phase difference between the right {I — +) and the left (/ = — ) layers of the same spin is always 
—Py as in a commensurate phase, while it is or Q'^ between up spin electrons and down spin electrons within 
the same layer. In other words the wavefunction proposed in Eq. (46) can only give a commensurate phase, because 
the effect of in-plane magnetic field has been automatically taken into account by transforming the layer index basis 
into the noninteracting energy eigenstate basis as shown in Section II B. More precisely, following Ref. [36], we can 
define three different states according to the phase difference of electrons in different layer and spin quantum states: 
(i) fully commensurate state, if Qx = Qx = 0? and hence spin up and spin down electrons in the same layer have the 
same winding phase determined by the in-plane field; (ii) partially commensurate/incommensurate state, if Qx 7^ or 
Q'^ ^ 0, and hence spin up and spin down electrons in the same layer have different winding frequency, although the 
phase difference between electrons of the same spin but in different layers still oscillates with a wavevector determined 

(1 2) (1 2^ 

by the in-plane field Py] (iii) fully incommensurate state, if Q_^Y{i) ~ ^-'lll) ^ ~^v^ ^^'^ therefore the tunneling 
energy becomes ineffective. According to Eq. (48), the fully incommensurate cannot be obtained from the trial 
wavefunction of Eq. (46), and only the fully commensurate and partially commensurate/incommensurate phases are 
the possible solutions. We note, however, that the fully incommensurate phase can be formally described by Eq. (46) 
if we use operators defined in Eq. (14) with Px = Py = and set As as = in the HF energy. Symmetry properties 
of the wavefunction, Eq. (46), may be easily described in analogy with Section III A3 (see Eqs. (43)-(45)). 



= 



2. Wavefunction for wide well systems 

For a wide well system at = 2N + 2, we have to separate the two kinds of level crossing possibilities in different 
regimes of the in-plane magnetic field: (i) the intersubband level crossing {W2) at small (see Figs. 1 and 2(c)), 
and (ii) the intrasubband level crossing (W2') at larger B|| (see Figs. 1 and 2(d)). In the first case, the intersubband 
level crossing for noninteracting electrons at even filling factors, there are two separate possible level crossings close to 
each other: ((1, 0), |) with ((0, N), t) and ((1,0), t) with ((0, N), l). The in-plane magnetic field (and hence Zeeman 
energy) is so small (for realistic situations in GaAs-based 2D systems) in this case that the two level crossings are 
actually very close, and therefore it is better to consider both of them simultaneously in a single trial wavefunction 
(constructed by the four degenerate levels) rather than consider them separately as two independent crossings. The 
simplest wavefunction in this case should be similar to Eq. (46), where the 4x4 matrix is block-diagonalized as two 
separate isospinors for each pair of the crossing levels. According to the symmetry arguments of the last section and 
the parity symmetry properties shown in Eq. (20) , such simple block-diagonalized 4x4 matrix representation of the 
trial wavefunction is further justified when considering a uniform phase only (i.e. Q± = Qj^ = 0, and Vfe and tp'j^ are 
constants), if the two crossing levels have different parity symmetries. As a result, we just consider the even-N case 
in this paper (so that level (1,0) and (0,N) have different parity and spin rotational symmetries), and speculate that 
the results for the odd-A^ case should be similar. Therefore, for the intersubband level crossing of a wide well system 
at even filling factors, we propose the following trial wavefunction, similar to Eq. (46): 
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^'l,+.fc 

l.-.k 



-.fc 



.jkQJl/2+i^/2 cos(^fe/2) e-*'='3»'o/2-i7/2 sin(^/>fc/2) 

.eifeQxi^/2+i7/2 sin(V'fe/2) e-''^'?- '0/2-^7/2 cos(V'fe/2) 










gikQ'jl/2+i'y' /2 cos(V'^/2) e-^'='5x'o/2-^772 sin(V';/2) 

_eifeQL'o/2+i772 sin(V';/2) e-^'=Q»'o/2-"'y'/2 cos(V'i^/2) 



'^ni,l,k-Qy/2 

n2,T,fe+Qy/2 
'^ni,T,fc-Q;/2 

'^n2,i,fe+Q;:,/2 

where ni = (1, 0) and ri2 = (0, iV) = (0, 2) are the Landau level indices of the crossing levels (see Fig. 2(c)). 

For the intrasubband level crossing at large -B|| region {W2'), only two noninteracting levels participate in level 
crossing (sec Figs. 2(d) and 5(a)): level 1 = ((0, 7V),|) and level 2 = {{0,N + l),t)j and all the other levels remain 
separated by a finite gap. In this situation one is allowed to consider only two degenerate levels when discussing the 
formation of a many-body state created by hybridization of levels 1 and 2. We will, however, still include mixing 
the next nearest levels, levels 3 = ((0, N), 1) and 4 = ((0, N + 1), |), in the theory, and consider the block-diagonal 
wavefunctions of the type given in Eq. (49) (but with ni = (0, A'') and n2 = (0, + 1)). In the lowest order in 
Coulomb interaction they agree with the two-level coupling wavefunction (c.f. Eq. (36)), but have an advantage 
that they allow us to discuss W2' level crossing point at the same footing as W2. The parity symmetry argument 
for a uniform phase can also be applied in this case, since levels (0, N) and (0, iV + 1) are always of different parity 
symmetry and the corresponding Zecman-split levels obviously have different spin polarizations. Therefore we believe 
Eq. (49) to be a reasonable trial wavefunction to study possible isospin winding and/or stripe order, although we 
cannot exclude the possibility that the ground state maybe stabilized by other more general wavefunctions. We do 
point out, however, that if we include the mixing of ((0, N), "f) with ((0, N + 1), J,), then ((0, A'' — 1), "f) mixing with 
((0, N + 2), I) arises exactly in the same order in Coulomb interaction, but with a numerically larger energy gap for 
the unperturbed levels. Therefore, unlike the double well system a.t 1^ = 4iV-|-2 (where the cyclotron resonance energy, 
co±, can be assumed to be much larger than the tunneling energy and the Zeeman splitting), considering only states 3 
and 4 and neglecting all the other levels in the wide well W2' case can not be energetically justified. In any case, our 
four-level coupling wavefunction agrees with the simple two-level coupling in the lowest order in Coulomb interaction, 
and we will show (in Section V) that its HF energy also gives the correct magnetoplasmon dispersion enabling us to 
study the possibility of a second order quantum phase transition to a state of broken spin symmetry [17]. Therefore 
we believe that Eq. (49) provides a reasonable trial wavefunction to describe the single wide well system for even 
filling factors — - the fact that the numerical implementation of the Hartree-Fock calculation using Eq. (49) is relative 
easy is an additional motivation to study it in details. 



(49) 



C. Stripe formation in the isospin coherent phcise 

In Section III A we have discussed some physical properties of a stripe phase, where tpk can be a periodic function 
of the guiding center coordinate, hence providing an oscillatory isospin polarization, (X^). A generic stripe phase 
discussed in this paper is provided by the hybridization of two crossing levels of different parity symmetries near the 
level crossing (or near degenerate) region, and therefore (at least) both parity (full space inversion) symmetry and 
translational symmetry are broken when a stripe phase is stabilized by Coulomb interaction. (Note that the spiral 
phase may also break translational symmetry without any stripe order; spin rotational symmetry around the total 
magnetic field direction may also be broken if the two coherent levels are of different spin directions in an even filling 
factor system.) Parity is broken in the non- fully polarized regions between the stripes that choose a spatial direction, 
and the translational symmetry is broken when the stripes choose their positions. Transitions between states of 
no-broken symmetries (i.e. fully isospin polarized states in our case) and states that break parity and translational 
symmetries simultaneously may happen in two ways. The first possibility is a direct first order transition when the in- 
plane magnetic field exceeds some critical value. To calculate the ground state wavefunction for this first order phase 
transition, we need to include many variational parameters in the theory (as shown in Fig. 6 and Appendix E) for the 
stripe phase function in addition to the spiral wavevector, Q±. The numerical calculation for this first order transition 
is very time-consuming in general. The second possibility is two consecutive transitions that break symmetries one 
by one: at the first transition (which could be either first order or second order) parity symmetry is broken through 
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a uniform superposition of the two crossing Landau levels, which have different parity symmetries, and no stripe 
order is present; at the second transition, the stripe order appears spontaneously with a concommitent breaking of 
the transitional symmetry via a second order phase transition. In this section we concentrate on the second scenario 
(two consecutive transitions) and develop a formalism for studying instabilities of an interlevel coherent phase toward 
the formation of stripe order. In some situations a spiral order may be stabilized in the first step, which, strictly 
speaking, breaks the translational symmetry even without a stripe order, see Section III A3. Appearance of the stripe 
order from such s spiral phase can also be described by the perturbation theory that we develop below as long as the 
stripe order appears via a second order phase transition. In Appendix D we will show that our perturbation method 
for probing the existence of a stripe phase is actually equivalent to studying the finite wavevcctor mode softening of 
a collective mode inside the uniform isospin coherent phase. We mention, however, that in general the perturbation 
calculation is easier to carry out than the mode softening calculation. 

The perturbation method consists of the following steps. First, we use the trial wavefunction in Eq. (36) to search 
for non-stripe phases (i.e. use Q*j_, 7, and ipk = V'o as the variational parameters) that minimize the energy of the 
system. If the optimal configuration has V'o = or tt, then no uniform many-body phases are stabilized near the level 
crossing (or near degeneracy) region. Formation of the stripe phase is still possible, but if it does occur, it happens 
via the first order transition, and we need to consider the explicit variational forms for the stripe phases, and compare 
their energies to the energies of the uniform isospin polarized phases (see Appendix E). Alternatively, if the optimal 
nonstripe configuration has ipo 0,Tr (and possibly finite Q*i_), the stripe phase may appear via a second order phase 
transition, which may be understood as the appearance of small oscillations in Therefore, we can choose the 
oscillation amplitude of tjjk to be an order parameter, and approximate ipk by the following formula: 



^fc = Vo +4Acos(A;g«^), 



(50) 



where |A| <c V'o 7^ Oj 9? the characteristic wavevector of the stripe (the stripe period a = 27r/g), is the only additional 
parameter we need in the perturbation theory. We then obtain the following expansion of Oj(g„) to the second order 
in A using the definition given in Eq. (37): 



ei(g„) = (5g„,o [em - 2A2cos(V'5)] - Asin(V'o*)K„,-q 
e2(g„) = (5,„,o [e^(0) + 2A2cos(V'S)] + Asin(V;o*)K„,-,^ 
Osiqn) = 5^^,o [em - 2A2sin(^*)] + A cos(Vo*)['5,„,-,^ 



A^ COs(V'o)[<5g„,-29 

■ A^ cos(V'o)[<5g„ -29 
A2 sin(V'o*)K„,-2q- 



^q^aq] 



-©(A^) 

-0(A3) 
■0(A3), 



(51) 
(52) 
(53) 



where 6*(0) {i = 1, 2, 3) are their extreme values. Putting Eqs. (51)-(53) in the expression of our HF energy as shown 
in the latter sections, we obtain the leading order (quadratic terms of A only) energy perturbation of a stripe phase 
from the HF energy of the uniform phase, E^mstripeii^o-: Q*l)- This result can be expressed as follows: 



(54) 



where the sign of EpJ^f{q; tpQ,Q*j_) determines the existence of a stripe phase: if the minimum value of Epj!^f{q; V'o; Ql) 
is negative and at a finite value of q = q*, we can claim that the original uniform phase is not energetically favorable 

compared to the stripe phase, and the groimd state can then be a stripe phase with iso(spin) winding vector Q*^ 
and stripe oscillation wavevector q* along x axis (we are using Landau gauge, A^yj{f)). On the other hand, if 

-^perf (9; V'o ' is positive for all q, then a stripe phase along y cannot be formed through a second order phase 
transition. If we want to study the possibility of stripe formation along x direction (i.e. stripe modulation is along y 
axis), we can do the same analysis as above, but using the Landau gauge, A[j.](r ). If both stripe phases are possible, 
we need to compare them and find the one of the lowest energy. 

Finally we note that the same approach can be also applied to study possible stripe formation via a continuous 
transition in the double well systems oX v = AN + 2 and wide well system at = 2 A'' + 2. In these cases we start with 
Eqs. (46) and (49) respectively, and use Eq. (50) and 



=Vo*+4A'cos(fcg%2)^ 
and expand the HF energy for small A and A': 

E^t^ipeiQ, q') = E^,l,,,,,M,,<*, Ql, Q'1.1 + [A, A'] . Ef,f,(g, q') 



A 
A' 



0{A% 



(55) 



(56) 



where E^fj(g, g') is a 2 x 2 matrix (we have suppressed all other fixed parameters, tpo, etc. for notational simplicity). 
Therefore if the lowest eigenvalue of E^^j(g, q') is negative and located at finite {q*,q'*), we can obtain a stripe phase 
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with total HF energy lower than the uniform coherent phases. On the other hand, if both eigenvalues of £^^^((7, q') 
are positive for the whole range of {<l,q') or its minimum value is at {q,q') = (0,0), then we conclude that no stripe 
phase should arise via a second order phase transition. Once again we emphasize that, in general, it is possible that 
stripe phases at large value of A are more favorable, for which the lowest order expansion in Eqs. (51)-(53) is not 
sufficient. But this would correspond to the first order transition to the stripe phase, and we do not have a better 
method to study its existence except for a direct numerical variational calculation (Appendix E). 



IV. WIDE WELL SYSTEMS AT z/ = 2iV + 1 



In this section we calculate the Hartree-Fock variational energy obtained by the trial wavefunction of Eq. (36), and 
show the numerical results for a wide well system at odd filling factors, v = 2N + 1. As mentioned earlier, there are 
two classes of level coherence in a wide well system (see Fig. 1): one is for the intersubband level crossing at smaller 
B|l (denoted by Wl), and the other one is the intrasubband "level near degeneracy" at larger B|| region (denoted 
by VFl'). By "level near degeneracy", wc mean a small (but non-zero) gap between energy levels of noninteracting 
electrons, which is much smaller than the Coulomb interaction energy. This is not a true level crossing (which would 
imply a zero gap rather than a "small" gap), but as we will show below, it is sufficient for inter level hybridization 
leading to non-trivial many-body ground states. For simplicity, we assume that the lowest 2N levels do not have any 
interlayer coherence while the spin polarized top level is allowed to have interlevel coherence with the lowest empty 
level of the same spin polarization. In Figs. 2(a) and (b) we show the corresponding quantum numbers for the relevant 
Landau levels of these two kinds of level coherence we consider in this section: in the small region, we consider 
the intersubband level crossing between ni = (1,0) and n2 = (0, iV), and in the large _B|| region, we consider the 
level near degeneracy between ni = (0, N) and n2 = (0, iV + 1). We note that this two level approximation is easily 
justified energetically in the Wl case, where a level crossing always ensures the noninteracting energy gap between 
the two crossing levels is smaller than than their energy separation with other levels. It is, however, less justifiable 
for the Wl' case, where the finite gap between the two coherent levels is just numerically smaller than their energy 
separation from other lower filled or higher empty levels. For simplicity, in this paper we will restrict our analysis in 
the Wl' case (intrasubband level near degeneracy) to A'' = (z^ = 1) only, and speculate that the results for other 
odd filling factors {N > 1) should be qualitatively similar. Note that in Fig. 1 there are more level crossings in the 
small region, e.g. crossing between (1,1) and (0,3), and also more crossing in the large -B|| region (not shown in 
the figure), e.g. the between (0, AT — 1) (spin down) and (0, A -I- 2) (spin up). For the sale of brevity, we will not 
discuss these additional level crossings in this paper. We believe the level crossings or level near degeneracy situations 
we consider in this paper are the most typical realistic ones for a wide well system, and the results for other level 
crossing situations should not be qualitatively different from the ones wc discuss in this paper. 

We use the trial state proposed in Eq. (36), and obtain the following expectation values: 

= ^si,S2'^fei,fe2^mi,m2 [cOS^(V'fei+Q^/2/2)(5^i,r-!i (5814/2 + sin^(V'fci-Q^/2/2)5mi,n2'5si,-l/2 + 5mifiO{N - m'l)] 
+^si,S2^si,l/25fei,fe2-Q,^mi,rJi^m2,n2e""^="''''"^''/^^'° sin(V'fe,_Q^/2/2) C0s(V'fe,_Q^/2/2) 6"'^ 

+^si,S2^si,l/25fe2,fei-Q,^mi,n2^m2,Sie''^'=^'''"'^''^^^'° ^H'4>k,-QJ2/'^) COsCV-fei -Q^/2/2) e'^ , (57) 

where fhi = {mi,m'^) (i = 1,2), and 6{x) is the Heaviside step function (= 1 if a; > and = otherwise). 



A. Hcirtree-Fock variational energy 

Using Eq. (57) we can calculate the single electron noninteracting energy from the noninteracting Hamiltonian of 
Eq. (6): 

Er = E^^m + E^e^iO) + J: 'E^, (58) 

where ^ ' means a summation over the core state. The Hartree (direct) energy per electron can also be obtained 
from the direct term of Eq. (8): 
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2n 



4E' [^nTn.,™™(0'0)ei(o) + fif^^,^^(o,o)e2(o)' 

Wm^u, (in, 0) e'^^] &s{qn - Q.)Qs{qn + Q.) 

{^^i'(9n,o)ei(g„)2 + £;^2i(g„,o)e2fe)' + 2i;^/(g„,o;Q^)ei(g„)e2((7„)} 

+2 [E^liO, O)ei(O) + E^iiO, 0)62(0)] + ^^^6 (9n, Q,)e3(g„ - Qx)' 
+'5q..o E ^nHqn, ; 27)e3((z„ - Q^JOsiqn + Qx) 

Qn 

+2^Q.,o E {^m ; 7)e3(g„ + Qx)ei{qn) + E^i ; 7)e3(g„ + Qx)e2(9„)} , 



e3(g„ + Qx)©2(gn)} 



(59) 



where in the last equation wo have introduced E^^^ [i = 1---6), and E^^^ (j = 1,2) to label the Hartree 

energies contributed by each corresponding term for the convenience of later discussion. Their definition is obvious 
from Eq. (59). Here g„ = 2wn/a is the stripe wavevector with a being the period of the stripe. Note that there 
arc three kinds of Hartree energies shown in Eq. (59): (i) the term E^^^ (i = 1 • • -G), which have no explicit phase 
(7) dependence, are finite for all value of (5± in general; (ii) the term £]^^, which has explicit e^^^ dependence, is 
nonzero only when Qy = 0; (iii) the terms E^j {j = 1,2), which have explicit e*''' dependence, are nonzero only when 
Qy =0 and for finite stripe order (when considering nonstripe phase, g„ = 0, V^^^ mm(^' ^) ~^)- Since our explicitly 
numerical results show that the third kind of contribution {E^j-) is always zero in the ground states we obtain near 
the level degeneracy region, we will neglect them throughout in our discussion. Comparing the E^-^ [i = 1 • • • 6) 
terms with the term, we find that their distinction arises from the fundamental difference of the ordered phases 
at finite and at = when no stripe order is present (note that when -ij^k is a constant, Qi{qn) = <5q„,o0i(O) 
and Qi{qn ± Qx) = '^g„,q=Q«0i(O) so that term is proportional to ^Q^fiSq^fi). For Q±_ finite, the state breaks 
translational symmetry and the invariance of energy with respect to 7 reflects a freedom of choice of the origin (see 
Eq. (44) and the discussion in Section III A3). It also signals the presence of a gapless Goldstone mode coming from 
the spontaneously broken continuous (translational) symmetry. On the other hand, for Q±^ = 0, the many-body state 
breaks only the discrete parity symmetry (see Eq. (43)). As a result there is an explicit dependence of energy on e*^''' 
with 7 = 0,7r being the two degenerate minima (see also the discussion below and in Section VIII D). The ground 
state selects either 7 = or tt via the Ising type transition, which describes the breaking of parity symmetry. No 
Goldstone mode exists in this case since the broken symmetry is discrete (i.e. Ising type). 
As for the Fock (exchange) energy per electron, we can have 

= ^ E ( ^nTm.mn. (9 ) E ^os{qnqyll)@i{qnf + Vl^,,n.n. {q) E cos{qr,qyll)&2{qnr 

q± I ^r. Qn 

(q ) E cos(g„(gy + Qy)ll)Qi{qn)Q2{qn) 

Qn 

+ 2E' [vZ^,rnnAq)Qm + VZn.,fnnA^)Q2{0) 

rn 

+2Vs^f^l,^^2f^2(9) cos{{q^Qy - qyQx)ll)^cos{qnqyll)Q3{qnf 

Qn 

+2SQ^,oRe [v^umM e''^] Y.cos{qnqyll)@3{qn - Qx)©3(g„ + Qx) 

Qn 
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Qn 

Pn ) 

= ^E{^^i'(9"'0)Qi(9")' + ^^2'(9n,0)e2(g„)' + 2£f3i(g„,0;Q^)ei(9„)e2(g„)} 

+ [i;f/(o,o)ei(o) + i;]?;i(o, 0)62(0)] +E^F6'('/n,0;O„Q,)e3(<?„)2 

+^Q.,o E ; 27)e3(9„ - Qx)e3(g„ + Qc.) 

+2'^Q.,o E {^Fiiln^ ; 7)03(«n + Qx)Qi{qn) + E%\qn, ; 7)e3(«n + Ox)e2(5„)} (60) 
gn 

where -E^^ (i = 1 • • • 6), f^^, and Ep^j^ {j = 1, 2) are introduced to label each contribution of the exchange energy. 
The dependence of these exchange terms on the phase 7 associated with possible broken symmetry behavior is the 
same as discussed earlier for the corresponding Hartrce energy terms. 

We think that it is worthwhile to emphasize again that Eqs. (59) and (60) are based on a specific choice of the Landau 
gauge for the vector potential, A^y]{f), in which the particle momentum is conserved along y axis (perpendicular to 
the in-plane magnetic field). To obtain the HF variational energy for a stripe phase along x direction, we can choose 
the alternate gauge, A^^j (r ) , in which particle momentum is conserved along x direction, to construct a many-body 
wavefunction similar to Eq. (36) (see details in Appendix A). If no stripe phase is stabilized, the results obtained 
in these two gauges arc identical. To save space, we will not show the HF variational energy obtained in the second 
gauge throughout this paper, although we take it into consideration in our numerical calculations. 

B. Magnetoplasmon excitations 

Before showing the results of minimizing the HF energy, it is instructive to address the close relationship between 
the HF variational energy shown in Eqs. (58)-(60) and the collective magnetoplasmon excitations of the conven- 
tional incompressible quantum Hall states (i.e. the isospin polarized states). In an integer quantum Hall system, 
magnetoplasmons are collective modes associated with magneto-exciton excitations above the Fermi energy that 
can be theoretically studied by using the generalized Hartree-Fock (or time-dependent Hartree-Fock) approximation 
[48,17,33], which is correct to the leading order of the ratio of the electron interaction energy to the noninteracting 
Landau energy separation. The softening of the magnetoplasmon mode indicates that the system may undergo a 
second order phase transition from a usual isospin polarized state (i.e. the uniform quantum Hall state) to a new 
symmetry-broken ground state, which is precisely the same as that obtained by minimizing the variational HF en- 
ergy shown in Eqs. (58)-(60). Moreover, the full analytical expression of the magnetoplasmon dispersion can be 
obtained from the uniform variational HF energy (i.e. same as Eqs. (58)- (60) but considering ipk = ipo or equiv- 
alently Oi{qn) = 0i(O)(Jq„,o) by taking small It/jqI expansion from the isospin up ground state or by taking small 
Itt — '00 1 expansion from the isospin down ground state. For example, if we consider the Wl case with isospin up state 
(ni = (1,0)) being the highest filled level, the HF energy of Eqs. (58)-(60) can be expanded to the leading order of 
82(0) (i.e. small Vo) to obtain (using 63(0)2 = 61(0)62(0)) the following total HF energy: 

E^ki^o) = Ell^ + E^^;l\Q^)Q,{Q) + 5Q.,o^Q„o^2;I,'(7)e2(0) + 0(62(0)2), (61) 

where the first term is the total electron energy of the isospin up state (the ground state), including the HF self-energy 

correction, the second term is the plasmon dispersion shown below, and the third one is the additional point energy 
shift associated with the broken parity symmetry. More explicitly we have 

En's = E°,'Z + \ [Ehi{^, 0) + £fii(0, 0)] + E^l{Q, 0) + iJf^^O, 0), (62) 

where the first term is the noninteracting energy, the second term is the self-energy produced by electrons within the 
top level (I is for double counting), and the third term is the self-energy produced by electrons in the core state (see 
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the definition of E^^p-^^{0,0) in Eqs. (59) and (60)). The magnetoplasmon excitation energy in the right hand side 
of Eq. (61) gives 

E'^TiQ^) = - + ^^!t - + <6 (Q- Qy) + 0; Q., Qy), (63) 

where E^^^ is the Hartree-Fock self-energy of level n and spin s. Eq. (63) is exactly the same as the magnetoplasmon 

excitation energy of the incompressible (isospin up) quantum Hall state obtained directly from the time-dependent 
Hartree-Fock approximation (TDHFA) [17]. The contributions from the bubble diagrams (the direct term) and from 
the ladder diagrams (the exchange term) correspond to the last two terms of Eq. (63) respectively. Note that 
the energy of the Q± = point is disconnected from the rest of the spectrum due to the last term in Eq. (61), 
Sp^'J^ij) = £j[^{0,0;-f) + f (0, ; 7) < 0. This reflects the fact that a many-body state at Q±^ = breaks only 

a discrete symmetry and should not have Goldstone modes. Alternatively continuous dispersion for Q± ^ would 
break the continuous translational symmetry and lead to Goldstone modes. The plasmon dispersion obtained are 
based on the uniform integer quantum Hall state, or equivalently, the isospin polarized ground state, therefore Eq. 
(63) has to be changed if we want to study the dispersion of the collective modes inside the symmetry-broken ground 
state. We mention that such a complete equivalence between the Hartree-Fock ground state energetic calculation and 
the corresponding collective mode dispersion follows from the Ward identities, and has also been used extensively in 
Ref. [33] in discussing the canted anti-ferromagnetic state in bilayer systems. 

C. Results I: intersubband level crossing {Wl case) 

In Fig. 7(a) we first show the energy dispersion of the magnetoplasmon mode (in charge channel only) obtained 
from Eq. (63) for u = 5, near the intersubband level crossing point at B*^ ~ 2.33 Tesla (after including the self-energy 

correction) for realistic GaAs 2D systems. We also use a filled circle to denote the excitation energy at q± = 0, that 
are disconnected from the rest of the spectrum by a negative energy shift, S^'^fi^, according to Eq. (63). In the system 
parameter range we consider here {B± = 3 T and the bare confinement energy, uo, is 7 meV), this disconnected energy 
shift is very small (< 0.01 meV). When the in-plane magnetic field is above B*^ = 2.25 Tesla, the magnetoplasmon 
mode is softened at qj_ = 0, indicating a second order phase transition toward a many-body coherent state breaking 
the parity symmetry. Strictly speaking, only the disconnected point at q± = is softened at and the whole 
collective mode dispersion will be modified for B|| > B*^ inside the new symmetry-broken phase. 

In Fig. 7(b), we show the HF energy calculated from Eqs. (58)-(60) aroimd the intersubband level crossing point. 
We find an isospin coherent phase (V'fc = V'o 7^ 0) with no spiral order (Q^ = 0) and no stripe order) in addition 
to the isospin polarized quantum Hall states within a small range of the in-plane magnetic field (2.25 < B|| < 2.40 
T) for the chosen system parameters. According to the symmetry analysis discussed in Section III A 3, the coherent 
phase only breaks the parity symmetry of the system and therefore has no Goldstone mode. This is consistent with 
the result studied by the mode softening of the collective excitations shown in Fig. 7(a). 

We note that the new coherent phase, breaking the discrete parity symmetry, is similar to the ferroelectric state 
observed in ferroelectric crystals [42]. More precisely, the electric dipole moment, (r ), is obviously zero if the ground 
state has a definite parity, while it can be nonzero if the ground state mixes two states of diff'erent parities. Therefore 
we think the simple coherent state we find above in the intersubband level crossing region of odd filling systems is a 
"ferroelectric" quantum Hall state with finite electric dipole moment. The recent experiments observing anomalies in 
the Shubnikov-de Haas oscillations of a wide parabolic well in the presence of a tilted magnetic field may be due to 
the existence of such coherent states [50], but more definite experimental work would be needed to settle this point. 

D. Results II: intrasubband level near degeneracy {Wl' case) 

In Fig. 8 we show a typical magnetoplasmon mode dispersion in the charge channel of a wide well at = 1 in 
the large region near the intrasubband level near degeneracy point, where the interaction energy is of the same 
order as the noninteracting energy separation. (Note that the bare confinement energy uo is 3 meV here.) When 
is larger than 30 Tesla, we find a mode softening at a finite wavevector perpendicular to the in-plane magnetic 
field (i.e. along y axis). In the same figure, we use filled squares, triangles, and circles to denote the energies of 
the zero momentum excitation, which is different from the long wavelength limit of the plasmon curve by an energy 
l^rd oi\o 1^*^)1 '^•^^ meV (see Eq. (61)). In sharp contrast to the Wl case shown in Fig. 7(a), the collective mode 
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softening occurs here at a finite wavevcctor, Q±^ = (0, iQp rather than at Q±^ = 0, showing an isospin spiral order 
in this system. Therefore the ground state can be a spiral phase if only one of the ordering wavcvectors, (0, ±Q*), is 
present, or it can be a coUinear spin density wave if there is an ordering at both wavevectors with equal amplitude. 
We have not been able to write a wavefunction for such a coUinear phase to compare its energy with the spiral phase, 
and therefore we cannot rule out the possibility that a coUinear phase can also be a ground state in the Wl' case. 

Now we have to investigate if such uniform coherent isospin spiral phase is stable against the formation of a stripe 
phase. We use the perturbation method developed in Section IIIC and calculate the perturbative energy £'^ff(g). In 
Fig. 9 we show our numerical results for EpJ^f.{q) as a function of q for several different values of By. Both gauges 
of the vector potential, A[x]{f) and A[y]{f), are considered in calculating Ep^^^{q) as indicated in the figure caption. 
When is larger than a critical value (it is also about 30 Tesla in this situation), the minimum of is located 

at a finite wavevector {q* ~ 0.3 x 10^ cm~^) along x axis, showing a stripe order with isospin 1^ modulating in the 
X direction with a period 2n/q* ~ 2000 A. Therefore, combining the two results above, we conclude that an isospin 
skyrmion stripe phase (sec Section III A) can be stabilized, with the stripe normal vector along x direction and the 
spiral winding vector along y direction. In our numerical calculation, we do not see signature for any intermediate 
phase (e.g. isospin spiral phase without stripe order) between the isospin polarized (incompressible) quantum Hall 
state and the isospin skyrmion stripe phase — the local minimum of Epert occurs at finite wavevcctor simultaneously 
with the plasmon mode softening. Therefore, following the results of Section III A 3, the spiral order breaks the 
translational symmetry along y direction, while the stripe order breaks the translational symmetry in x direction 
(parallel to the in-plane field). As discussed in our earlier paper [8], such skyrmion stripe has finite topological isospin 
density that leads to charge stripe order with stripes perpendicular to the in-plane field. This should lead to anisotropy 
in charge transport with larger conductivity along the stripes, i.e. perpendicular to -B|| . 

In Fig. 10 we show the phase diagram of the wide well system at = 1 in a strong in-plane field. The usual 
incompressible integer quantum Hall state is favored at small well width (large bare confinement energy) and/or small 
values. At larger well width and/or stronger _B|j field, the system undergoes a second order phase transition toward 
an isospin skyrmion stripe phase with translational symmetries broken in both x and y directions (parity symmetry 
is of course also broken). In extremely large By and large well width, we expect the isospin skyrmion stripe phase to 
evolve toward the Wigner crystal phase, which, however, is not included in our present theory. 



V. SINGLE WIDE WELL SYSTEM AT = 2Ar 2 



For a wide well system at even filling factor, ~ 2N + 2, we also consider two kinds of level crossings (see Figs. 
1, 2(c) and (d)) : one is the intersubband level crossing between levels ((1,0), J, (t)) and ((0,7V),t (i)) in the small 
By region (14^2), and the other is the intrasubband level crossing between levels ((0, N), l) and ((0, N + 1), t)) in the 
large S|| region {W2'). For simplicity, we do not discuss the system behavior at yet higher fields, for example, when 
there is a crossing between levels ( (0 , A'' — 1 ) , J, ) and {{0,N + 2),^). The main difference between a level crossing in an 
odd filing system and the one in an even filling system is the spin degree of freedom. In the odd filling situation, the 
two crossing levels are of the same spin so that only isospin degree of freedom affects the existence of a novel coherent 
phase via their different orbital wavefunctions (and the spin degree of freedom is essentially frozen). In even filling 
situation, however, the two crossing levels are of opposite spin polarization, so that a stabilized many-body coherent 
phase must also break the spin rotational symmetry. This fact leads important consequences, since spin and isospin 
arc not equivalent in their roles: the Coulomb interaction docs not flip spin polarization but may flip the isospin 
polarization of each scattered electron (see Fig. 5(b)). We will discuss this subject in more details in Section VIII B. 

As shown in Section IIIB2, when discussing possible many-body states around level crossings at v = 2N + 2, 
we will consider trial wavefunctions that mix the four closest levels around the Fermi energy, the two highest filled 
Landau levels and two lowest empty Landau levels, and will assume that the lower 2N (core) levels are completely 
filled (frozen) and do not participate in the coherent hybridization process. Using Eq. (49), we obtain the following 
expectation value similar to Eq. (57): 

= <5mi,m2'5cri,<72 4i,fe2 [cOs'^{lpk+Qy/2/2)5rniMi6cri,-l/2 + sin^ (V'fc-Qy /2 /2)/2)(5^i ,^2 (5^^ ,1/2 

+ COS^(^fe+Q;/2/2)<5mi,ni 1^^1,1/2 + sin^(V'fc-Q^/2/2)<5mi,n2'^<Ti -1/2 + SmiflOiN - m[) 
+^mi,n,'^rn2,n2'5^:,-l/2'5a2,l/24i,fc2-Q,e-^Q==('=^-«''/2)'o sin(V'fe,_Q^/2/2) C0s(V'fc2 -Q„/2/2) C'^^/^ 
+^mi,n2'5m2,ni'5ai,l/2<5a2,-l/2<^fc2,fei-Q«e*'3-('=^-^^/2)'o sin(V;fei _Q^/2/2) C0s(V'fei -Q„/2/2) 6^^/^ 
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+<5A„n.^™„n,5.„-i/25.„i/25fe,fe,-Q;e^«i('=-«;/2)'o sm{i,'^^_Q,,j2) cos(V;;,_Q, /2/2) e^^'/2^ (64) 
where ni = (1, 0) and n2 = (0, TV) for W2 case, while ni = (0, N) and ^2 = (0, TV + 1) for W2' case. 

A. Hcirtree-Fock variational energy 

The noninteracting single electron energy can be obtained from Eq. (11): 

Er = El'Z&m + eIZQ2{0) + El'^,e[{0) + EZe'M + ^ 'Eff, (65) 

where the last term is a constant energy shift from the frozen core states. The Hartree (direct) and the Fock (exchange) 
energies per electron are respectively 

= ^ E {v^nMAQn) [eiiqnf + e'M'] + V^n.,n.M [e2{qnr + ©2(9n)'] 

+2V^nA.«ini(9«)ei(g„)e'i((;„)cos(5„(g, - o;)/g/2) 

+2Varn.,n.n.ten)e2(5„)e'2(9„)cOs(9„(Q, - Q'y)ll/2) 

+'^yZn,,n.nMn) [ei(9„)e^(g„) + e2(g„)e;(g„)] cos(9„(Q^ + Q'y)ll/2) 



Vl,n,,inrn{^) (©1(0) + 61(0)) + W,n.,inrn{^) (©2(0) + 6^(0)) + ^ 'V^^,^,^,{Q) 



(66) 



and 

+'^V^n,,n,nA?±) E COs(9n5,/g) COs{qn{Qy + Q'y)ll/2) [BM&M + e2{qn)e[{qn)] 



+ 2E'^nT™(9l) [e»i(o) + 0'i(o)] + [02(0) + e^(o)] + 2 ^ FX^„^,^,(g-) 

m m mi,m2 

_2 

+^ E ^nTsi,S2n2 cos{{Q^qy - Qyqx)ll) ^ cos(g„g;,/g)e3(5„)2 
_2 _ , 

+ E ^«Tni,n2n2(9J-) COs((Q;,g;, - Q;,gr.)^o) E '=°^(^"«2''o)03(9n)^ 



2f2L 



E^e&."i«.(«^) e^^'^^'^'^] E03(9n - QMiln + Q'Jcos{{q,Qy - qnqy)ll), (67) 



2f^_L 

where the last term in the exchange energy shows the interplay between the two isospinors as mentioned in Section 

TUB 1, and this term is nonzero only when Q'^ = ~Qi_ (for simplicity, wc have chosen their stripe wavcvcctors, g„ 
and q'^, to be the same). Note that this is also the only term which depends on the phase, 7 + 7', in the Hartree- Fock 
energy. This is because, in the even filling systems, the two crossing levels are of different spin polarizations. Therefore, 
spin symmetry breaking of the coherent state gives a continuous energy degeneracy (i.e. 7 — 7' is arbitrary), while 
the breaking of parity symmetry selects 7 + 7' = 2m7r, where m is an integer. 
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B. Magnetoplasmon excitations 



Using the arguments similar to those in Section IV B, we can also obtain the magnetoplasmon excitation energy 

of the oven filling system in the usual incompressible quantum Hall ground state by taking small ©2(0) and ©2(0) 
limits in the above Hartree-Fock variational energy, Eqs. (65)-(67). The result is equivalent to solving the eigenvalue 
problem of the following 2x2 matrix: 



(Q±) 



n2T)"ii 



n2T."ii 



(Q±) 



T^W2 

Ex' 



"21, ml 



n2i,niT 



Ex 



iQ±) 



(68) 



jpO,W 
n2,T(i) 



n2,T(i) 



where we have used Q'. = -Qj_; AE°'Z,^ - ut^ = E°'Z,:^ - ^^°'^rt^ and AE?-^;!^, ,,,, 
"^fh^d) ~ ^'"^ respectively the noninteracting energy and the HF self-energy difference between the two 

relevant levels. The definition and the explicit expression of the HF self-energies, E^^^'^^, are similar to those in the 
odd filling systems. The two electron- hole (exciton) binding energies are respectively 



(69) 



and 



E^:\Ql_) = ^E^nTn2,mn2(9^)cOs((5.Q,-«,Qx)/^). 



(70) 



We note that Eq. (68) is exactly the same as the magnetoplasmon dispersion matrix derived for the triplet spin 
channel in the time-dependent-Hartree-Fock approximation [17], demonstrating that our proposed four level trial 
wavefunction. Eq. (49), is adequate in investigating the existence of new broken symmetry phases. On the other 
hand, we note that the other two spin singlet magnetoplasmon modes [17] cannot be obtained in our theory, since the 
trial wavefunction of Eq. (49) is still not of the most general form for the four level degeneracy. As mentioned earlier, 
this fact will not affect any of our results or conclusions shown in this section, because these two singlet excitations 
are relatively higher energy excitations and are of different symmetries from the lowest one we consider here. For the 
purpose of understanding quantum phase transitions in the system, it is crucial to have the correct description for the 
low energy sector of the relevant Hilbert space, and clearly our four level trial wavefunction of Eq. (49) accomplished 
that very well. Since the 2D magnetoplasmon dispersion calculation of the integer quantum Hall system has been 
reported before [17], we will not further discuss the magnetoplasmon dispersion and just focus on the HF variational 
energy calculation. 



C. Results I: intersubband level crossing {W2 case) 



For level crossing in the small B\\ region, our numerical calculation shows that there is no many-body coherent 
phase with total energy lower than the (uniform) isospin polarized states within our Hartree-Fock approximation. In 
other words, we find that such level crossing always introduces a (trivial) first order phase transition with a sharp 
polarization change in the narrow region of level crossing tuned by B\\ (see Fig. 1). This is related to the resistance 
hysteresis recently observed and discussed in Refs. [51,52] — small domain walls may occur during the first order 
phase transition separating the two polarizations so that the resistance shows a hysteretic behavior when the external 
electric field is swapt. Although such domain wall physics associated with intersubband level crossing induced first 
order transition is of intrinsic interest, we do not include this possibility in our theory since our interest here is to 
classify the (second order) quantum phase transitions between nontrivial quantum Hall phases (see a brief discussion 
in Section VIII E. 



D. Results II: intrasubband level crossing {W2' case) 

For the level crossing at large B\\ region, only one level crossing between Hi = (0, N) of spin down and n2 = (0, A'^+1) 
of spin up occurs and the next nearest two levels do not cross in the noninteracting energy spectrum. From Eqs. 
(65)-(67) we find that if we fix tl)). = tjjo to be uniform and finite {0 < tjjk < tt), the HF energy is always minimized 
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at a finite winding wavevcctor along y direction, Q± = ±(0,Q*) with Q* ~ 0.75 ?g ^ (the magnetic field is along x 
axis), so the optimum state of broken spin symmetry must have a spiral or a coUinear spin order (the latter happens 
when ±Q components are present simultaneously). However, in optimizing the HF energy with respect to tpo we 
find that the minimum is always at iJjq = or tt, so states with broken spin symmetry are not favored at the HF 
level. In Fig. If we compare the HF energies of the isospin spiral, isospin spiral stripe, and isospin skyrmion stripe 
phases, calculated using typical realistic system parameters (for GaAs 2D systems) and by employing a more general 
phase function, tpk, as shown in Fig. 6. For the sake of comparison we fix most parameters in the phase function of 
tpk and let 4> = tpi — 1^2 to be the only free parameter for the stripe phases (we have tried different ranges of these 
variational parameters, but the results are qualitatively similar to Fig. 11 and no exotic many-body phase is found). 
In the horizontal axis of Fig. 11, we define the spin polarization to be ©2(0), which is proportional to the density 
of spin triplet excitons excited from the top filled level to the lowest empty level. Another coherence parameter, 
©2(0), is chosen to be zero in Fig. 11 because the results do not depend qualitatively on this choice. In Fig. 11 
we choose = 11 T in the calculation, slightly lower than the level crossing point at S| = 11.1 Tesla (note that 
this is the renormalized level crossing field including the HF self-energy correction and is therefore lower than the 
noninteracting result). We find that, within our Hartrec-Fock calculation, the energy of the spiral phase is a convex 
curve as a function of the polarization, and therefore has no energy minimum between ©2(0) = (spin unpolarized 
state) and ©2(0) = 1 (spin polarized state), i.e. either the spin unpolarized or the spin polarized state is always 
the lowest energy ground state, depending on whether the magnetic field is smaller or larger than B*^ = 11.1 Tesla. 

Note that we have chosen the spiral wavevector to be = (0,Q*), which is the extreme value for the lowest HF 
energy for < ©2(0) < 1. Unsurprisingly, we find that this Q* is the same as the wavevector obtained from the 
near softening point of the magnetoplasmon excitations [17]. (Therefore; we can exclude the simple (incommensurate) 
isospin coherent (stripe) phases from the energy comparison, because they do not have spiral order and must have 
higher energy than the three cases shown in Fig. 11).) Among the three many-body states of spiral order, we find 
that the spiral phase always has the lowest energy. Our HF calculation therefore suggests a first order transition 
between spin polarized and unpolarized ground states with no intermediate broken symmetry phases in between (at 
least within the HF theory). We find, however, that the HF energy differences between simple incompressible states 
and exotic many-body states, such as the spiral phase and the skyrmion stripe phase, are very small (< 0.1 meV). We 
therefore suggest (and speculate) that effects not included in our analysis, e.g. self-consistent calculation for the single 
electron wavefunction, lower level screening, and/or nonparabolic effects of the well confinement etc., may now well 
stabilize the many-body states. It is also plausible, given the smallness in the HF energy difference, that higher-order 
corrections beyond the HF theory could stabilize exotic quantum order in this situation. 

In Ref. [15], Pan et. al. observed strong anisotropic longitudinal resistance when the in-plane field exceeded a 
certain critical value. This is very suggestive of the skyrmion stripe phase, since the latter has charge modulation in 
addition to the spin modulation, and therefore should lead to strong transport anisotropy. It is useful to point out 
that the direction of the charge modulation in the skyrmion stripe phase is fixed by the applied parallel magnetic 
field: the winding of the transverse components of spin, i^x,y)i is set by Q±_ and is perpendicular to B\\, while the 
modulation of (X^) is along By, so we have effectively a one-dimensional charge density wave that goes along By (see 
also the discussion about skyrmion stripe phase in Section III A 2) . Hence the expected "low" resistance direction of 
the skyrmion stripe phase is perpendicular to B\\, which is what was observed in Ref. [15]. 

An alternative interpretation of the resistance anisotropy has been recently suggested in Ref. [53], where Chalker 
et. al. argued that surface disorder will form domains close to the first order phase transition, that have anisotropic 
shape due to the presence of the tilted magnetic field. Naive argument would suggest that these domains differ only in 
the spin structure and should not contribute appreciably to the transport anisotropy. Analysis presented in this paper 
suggests, however, that boundaries between different domains in this case should be accompanied by the topological 
spin density, which leads to change density modulation and may lead to large resistance anisotropy. In addition there 
is always the possibility (already mentioned in Section V C) that a direct first order phase transition from the spin 
unpolarized to the spin polarized phase will give rise to domains with different (up or down) spin polarizations in the 
(effectively Ising) ferromagnetic phase. Again, these domains, separated by domain walls, would differ only in the spin 
orientations, and it is unclear how this could give rise to the observed resistance anisotropy seen in the experiments 
[15,16]. Also, the domain structure should lead to hysteretic behavior in the observed resistance, which has not been 
reported. We emphasize, however, that the possibility of a direct first order (Ising type) transition in the experiment 
of Ref. [15] cannot be ruled out — in fact, our HF calculation does indeed predict such a transition (but with very 
fragile energetics) as discussed above. 

To summarize, within our approximations we do not find any exotic phases as a true ground state close to the 
level crossing point of 1^2'. However, the exotic phases are very close in energy, and therefore we speculate that the 
resistance anisotropy observed in Ref. [15] may arise from the skyrmion stripe phase, which could be stabilized by 
effects not included in our theory. 
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VI. DOUBLE WELL SYSTEM AT v = AN + l {Dl CASE) 



As mentioned in the Introduction, in the double well system at odd filling factor i/ = AN + 1, many interesting 

phenomena have been explored, such as the intcrlaycr coherence (and the commensuratc-incommensiiratc phase 
transition) [22,26,54,31], unidirectional charge density wave (stripe) state [55,24,46,56], and the in-plane magnetic 
field induced charge imbalance phase [57] of the odd filling systems. Following our earlier work in Ref. [8], we will 
use a general trial wavefunction, Eq. (36), to include the isospin stripe order and the spiral order simultaneously, to 
obtain a rich quantum phase diagram within a single unified theory including all the effects mentioned above (which 
in the past have been studied in separate works using different techniques). In this section, we assiimc the following 
(reasonable) ordering of energy scales for the double well system: ui± ^ lo^ ^ ^sas (see Eq. (15)), so that the lowest 
AN filled levels can be assumed to be frozen core states with no coherence effect, and only the top filled level has 
coherence with the lowest empty level of the same spin polarization (but opposite parity) . Within this approximation 
it is more convenient to follow the standard convention in the literature and consider the isospin within the layer 
index basis to calculate the HF energy. 

For the convenience of later discussion of the symmetry properties of the many-body wavefunction, we first write 
down the many-body state explicitly in the layer index {I = ±) basis (for notational simplicity, we suppress the 
Landau level index and the spin index throughout this section): 

\^E\i^>^,Q±,l)) = U{^'^'''''°^'^^"^osi^^^^^^^ (71) 

k 

We then obtain the following expectation value for the HF energy calculation: 

= Si.^iJkiM [cos^{'^ki+Qy/2/2)Si^^+i+sin^{ipki-Qy/2/2)/2)Si,-i] 

+Si,,.i, [cos(Vfc,-Q,/2/2)sin(Vfc,_Q^/2/2)4,+Q,/2,fc.-Q,/2'^i„+ie-"5-('=^-««/')'°-*'' 

+ sin(Vfa-Q,/2/2)cos(Vfa_Q,/2/2)4.-Q,/2,fo+Q,/2'5i„-ie''3»^'=^-«''/'^'«+'^l. (72) 



A. Hartree-Fock variational energy 



Neglecting the constant energies associated with Landau levels and Zeeman splittings, the noninteracting single 
electron energy is entirely the tunneling energy: 



^0 



= ~ XI ^ \e~^''^y^oe-''^-''^°-'^ + e''=^«'oe''3-'='o+»7 Sp^^g^ sin(V'fe/2) cosiipk/^) 



= -AsAsSp^,Q^Q3iPy + Qx) COs{-f). 

The Hartree and the Fock energies can be written as (obtained by using Eq. (13)): 

Ns 



(73) 



= ^ E \^NN,NN{Qn) {Ql{qnf + e2{qnf) + 2V^}^^^M C0s{qnQyll)eMQ2{qn)] , (74) 



and 



= E ^ E {^Nk^,NNiQ) COsiqr^qyll) [@l{qnf + Q2{qny 



+ '^^NNNwiq) cosiiQxqy - QvQx)li) cos(g 

nqylo)^3{qn) 



where we have used V^I^Jfj^{qn) = VNN,NNiqn), and neglected the constant energies associated with the frozen core 
levels. Combining Eqs. (73)-(75) we obtain the total HF energy as follows: 



E^j, = -AsAsSp^,Qye3{Py + Qx) cos(7) + ^ E { [^Uqn) + E+{q„)] ■ {Qiiqnf + e2{qnf) 
+2E],{qn;Qy)ei{qn)Q2{qn) + 2Ep{qn;Q±)@3{qnf] , 



(76) 



where the definition is obvious by comparing Eq. (76) with Eqs. (73)- (75). The only 7-dcpendent term is from the 
tunneling amplitude, reflecting the fact that the isospin rotational symmetry is broken by electron tunneling. We will 
discuss the symmetry properties in details later. 
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B. Commensurate, incommensurate, and chcirge imbalance phases 

We first analytically discuss a special class of many-body states implied by Eq. (71) in the absence of any stripe 
order, i.e. Vfe = V'o ^ constant. Taking = in Eq. (76), we obtain the following simplified HF energy (7 is set 
to be zero): 

E^^'iQi.) = -63(0) [AsasSq,,pJq^,-p, + Ea{Q^)Qs{0)} , (77) 

where Ea{Q±) = E+{0) - E^{0]Qy) + £^+(0) - Ep{Q;Qi_). We have used 61(0)62(0) = 63(0)^ and neglected 
the irrelevant constant energy. Following the existing literature, we separate the discussion in two parts: systems in 
an incommensurate state (i.e. the coherent phase we defined in Section 111 A 2 with the extreme value of winding 
wavevector, Q\ = 0), and systems in a commensurate state (i.e. the spiral phase defined in Section III A 2 with 
= {—Py,Px))- When the system is in an incommensurate state, the tunneling amplitude is effectively zero 
according to Eq. (77), and therefore the extreme value of -(/iq is determined by the sign of Ea{Q*i)- For Ea{Q*i) > 
(i.e. the Hartree energy dominates the Fock energy), the minimum energy is at 63(0) = 1/2 or = w/2, indicating 
an equal population of electrons in the two layers. However, if Ea{Q*j_) < (i.e. the Fock energy dominates the 
Hartree energy), the minimum value of E^^iQl) is at 6^(0) = (i.e. Vo = or tt), and therefore we obtain a 
fully spontaneous charge imbalanced state [57], where all electrons like to accumulate in a single layer rather than 
distribute equally in the two layers (this is true, of course, only within the HF approximation where the correlation 
energy is totally neglected). In the commensurate state (i.e. Qx = —Py and Qy = P^) and £Ja(Q1) < —^sas, the 
total energy in Eq. (77) is minimized at < 63(0) = \ Asas /"^EaIQ*^)] < 1/2, showing a spontaneous partial charge 
imbalance phase. Otherwise, for Ea{Q±) > —Asas, the commensurate phase always has equal number of electrons 
in the two layers, i.e. the usual isospin paramagnetic phase. Therefore, the in-plane magnetic field can cause not 
only a commensurate-incommensurate phase transition, but also a spontaneous charge imbalance phase when the 
exchange energy is large. We note that such exchange driven spontaneous charge imbalance phases could arise even 
in the zero-field (i.e. Btot = 0) non-quantum-Hall bilayer 2D systems within a restricted HF approximation [68], but 
in the zero-field case the corresponding XY "isospin magnetic" state has been shown to be lower in the energy than 
the Ising-type charge imbalance phase for the long range Coulomb interaction. 

In the present paper, we only consider the long-range Coulomb interaction, which gives E^{0) — E^{Q\ Qy) = e^/e^g 
(e is the dielectric constant of the system), and therefore the Hartree electrostatic energy always dominates the Fock 
exchange energy (i.e. -Ea(Q_l) > 0), eliminating any spontaneous charge imbalance between the two layers. On the 
other hand, as will be shown later, we may obtain a commensurate stripe phase with a longer period (a 3> ^o) in the 
small layer separation region, which is the asymptotic behavior of the charge imbalance phase recently discussed by 
Radzihovsky et. al. [57]. 

C. Numerical results and the stripe phases 

Fig. 12(a) and (b) are the phase diagrams we obtain at zero temperature for i/ = 5 for the in-plane magnetic field 
fixed in x direction. Using the isospin many-body phases defined in Section III, phase I is the coherent phase; phase 
II is the spiral phase, where the optimal spiral winding wavevector Q*^ = (0,(5*) = {0,Px) is perpendicular to the 
in-plane field direction; phase HI is the coherent stripe phase, where the stripe direction can be in arbitrary direction, 
and phases IV and V correspond to the spiral stripe phase, where Q*j_ is perpendicular to B|| and the stripe is aligned 

in X direction (i.e. (Tz) modulates in y direction, parallel to Q). To obtain the stripe phase (phases HI, IV and 
V), we numerically minimize the HF variational energy by using a general stripe phase function as shown in Fig. 6 
and in Appendix E. (The perturbation method developed in Section HI C also gives similar results, since the stripe 
formation in this system is a second order transition.) As mentioned in the beginning of Section III, to get the energy 
of a spiral stripe phase with stripe normal direction perpendicular to the in-plane field, we can just rotate the in-plane 
field direction from the x to the y axis. Therefore the results shown in 12(a) and (b) are gauge independent. Note 
that the stripe period of phase V is very large (a ^ Iq), showing an asymptotic behavior of charge imbalance phase 
modified by the long-range Coulomb interaction. Using the existing terminology of the literature [46,22,26], phases I 
and II are the incommensurate and commensurate phases respectively, and phases HI and IV are the incommensurate 
and commensurate stripe phases respectively. We will use this terminology as well as the isospin phases mentioned 
above (defined in Section HI A 2) for later discussion in this paper. 
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We note, however, that our HF calculation does not incorporate the possibility of two decoupled compressible 
V = 2N + 1/2 states in each layer, which could be energetically favored for smaller N (lower values of v, e.e. ly — 1). 
We expect that the phases discussed in this paper arc more likely to be found for N > 0, when each of the layers 
becomes susceptible to forming a stripe phase [1,4,5], especially in the presence of a parallel magnetic field. Particle- 
hole symmetry implies that similar states should also occur at filling factors 4iV + 3 by interchanging the role of holes 
and electrons in the top filled Landau level. 



D. Symmetry properties of the commensurate and incommensurate states 

Applying Eqs. (21) and (31)-(32) to the wavefunction, <91, 7*)) in Eq. (71) (where i'l,Q*_i, and 7* denote 

the extreme values to minimize the HF energy), we obtain 

fc " " 

= i-if'^l^G^r-k + Ql, -7* + ^)> (78) 
^x(i?.)|*G(^fc,Ql,7*)) = n (e'^''^^'°+"*-''^''"^/'cos(V^/2)c^;t^^^/^,_Q./, 

k 

^,-K^Q:il^.-n.PJ2)n sin(^^/2)<'t^^^/,.^^./ J ILL) 

= \'^G{rk-R.Jil,Ql.l* ~R.{Ql + Py))). (79) 
%{Ry)\^G{rkMl,ll) = n (e^('=«^'o+7+«.P.)/2g.(fc-Q;/2)ii„ ^os{rj2)c1'l_Q,f, 

k " 

^Q-i{kQ'Jl+i+RyP./2)/2^i(k+Ql/2)R, sin(V;^/2)cf j^Q./^) |LL) 

= \^G{rk,Qln* -Ry{Ql-P.))), (80) 

where we have used the fact that the guiding center coordinates, k, can be shifted and Ylik ^ = in a completely filled 
Landau level. 

Now wc can study the symmetry properties of the many-body phases obtained earlier in the HF approximation. 
In-plane field, B\\, is fixed along x axis and hence Py = 0. According to Eqs. (78)-(80), we find that (i) for an 
incommensurate state (= coherent phase with neither isospin spiral nor stripe order: xp^ = 7r/2, Q\ = 0, and 7* = 
arbitrary), V\^q^) = T'xl'i'Q^) = l^'^"'^) and therefore the parity symmetry and the translational symmetry in x 

direction arc not broken. However, since 7J,(i?y)|^'G(^^, Q^, 7*)) = |4'G(V'fei Ql) 7* + RyPx)) is in general not equal 
to the original wavefunction, and the HF energy is independent of phase 7* (the tunneling term is effective absence in 
the incommensurate phase) , the translational symmetry in y direction is therefore broken, with an oscillation period 
of2n/\Px\. Note that the above translational symmetry breaking in the incommensurate phase is obtained within the 
mean field HF approximation, which is not the same as the soliton lattice ground state obtained in Ref. [31] by using 
an effective field theory to go beyond the mean field approximation. Wc should expect to see a gapless Goldstonc 
mode associated with the broken translational symmetry in the time-dependent Hartree-Fock approximation, (ii) 
For a commensurate phase (= spiral phase: ipl = 7r/2, Q± = {0,Px), and 7* = 0), we find that no symmetry is 
broken at all, since no distinct wavefunction is obtained after applying P, and Ty operators. In other words, the 
commensurate (or spiral) phase is basically the same as an isospin polarized state from the symmetry point of view. 
It is not surprising because one can show that such commensurate (or isospin spiral) state is essentially the same as 
the noninteracting ground state, Hfe ^li + s k\-^^) (where a]^ ^ ^ j, is defined in Eq. (14) for the noninteracting energy 
eigenstate), of the double well system, which is certainly an eigenstate of all of these symmetry operators, (iii) For 
an incommensurate stripe phase (= coherent stripe: ipl. = 4'k^a/2-Ki^^ Q-*- ~ ^' ^^'^ ~ arbitrary), we find that the 
translational symmetries in both x and y directions are broken. However, we note that another incommensurate stripe 
phase, whose isospin z component, (X^), modulates in y direction (perpendicular to i3||), is also a degenerate ground 
state at the mean field level, breaking the translational symmetry in y direction only. The wavefunction of this second 
incommensurate stripe phase cannot be simply described in the present Landau gauge, A[y] (r ), and therefore its energy 
is not shown in the equations presented in this section. We may, however, consider an equivalent state by assuming 
the in-plane magnetic field to be along y direction with stripe modulation along x so that the translational symmetry 
is broken only in the direction perpendicular to the in-plane field direction in this second type of incommensurate 
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stripe phase. Parity symmetry is, however, broken in both the incommensurate stripe phases, (iv) Finally, for a 
commensurate stripe (= spiral stripe: V^, modulates in y direction, Q±_ = {Q,Px), and 7* = 0), its wavefunction 
cannot simply be described by the present gauge A[j,](f), and we can obtain it by effectively rotating the in-plane 
field direction as described above. By changing the in-plane field direction from x to ij, the spiral wavevector becomes 
(5j_ = {—Py.Q), and Eqs. (78)-(80) tell us that both parity symmetry and translational symmetry in x direction 
(now it is perpendicular to the in-plane field) are broken, while translational symmetry in y direction (parallel to 
the in-plane field) is preserved. We then conclude that the parity symmetry and the translational symmetry in the 
direction perpendicular to the in-plane field (i.e. in y axis) are both broken in this commensurate stripe phase in the 
presence of field (in x direction). This direction of broken translational symmetry is also parallel to the direction 
of isospin spiral wavefunction and the stripe oscillation (as defined for the spiral stripe phase in Section III A3). 



VII. DOUBLE WELL SYSTEM AT v = m + 2 {D2 CASE) 



For a double well system at even filling factor, v = AN + 2, both spin and layer indices are involved in the level 
crossing region (see Fig. 2(f)), where the cyclotron resonance energy, u)±, is (realistically) assumed to be much larger 
than the tunneling energy {Asas) and the Zeeman energy (cj^) so that the two top filled levels and the two lowest 
empty levels belong to the same orbital quantum number, well-separated from all other filled or empty levels, which 
are considered frozen and neglected from our consideration. As mentioned in Section II B, an appropriate basis for 
the isospinor for this system is the basis of noninteracting energy eigenstatcs, which have definite parity and spin 
symmetries. (Our definition of isospin coherent and isospin spiral phases in Section III A 2 then just corresponds to the 
fully commensurate and partially incommensurate phases in the literature [36] . For the convenience of comparison, 
we will use the conventional language for most of our discussion in this section.) We note that several interesting 
phenomena associated with the CAF phase in this system have been studied in the literature in the recent years 
[33-39]. A spin symmetry broken canted antiferromagnetic phase can be stabilized in this system in addition to the 
usual symmetric and fully ferromagnetic states in the = situation [33]. In the canted phase electrons in the two 
layers hold the same z component of spin polarization, while they have opposite spin direction in the x — y plane, 
showing a two-dimensional antiferromagnetic order. When an in-plane magnetic field is applied, such a canted phase 
remaining in a commensurate state and docs not exhibit a commensurate-incommensurate phase transition at large 
B|| as observed in v = AN + 1 case [36]. However, the possible stripe phase formation in the presence of an in-plane 
magnetic field for a larger layer separation has not yet been explored in the literature. For the sake of completeness, 
we will use the trial wavefunction proposed in Section III B and the perturbation method developed in Section III C 
to study the possibility of a stripe formation in the double well system at even filling factors. 

For a double well system at v = AN + 2, we have argued that the trial wavefunction proposed in Eq. (46), although 
not the most possible wavefunction, should be a reasonable approximation to describe the ground state in the presence 
of in-plane magnetic field. Thus we again start from calculating the following expectation value by using the trial 
wavefunction, shown in Eq. (46): 

(*G^ ,ai ,0-1 ,fci '^ma ,Q2 ,0-2 ,fc2 I *G^) 

k2 [cos {lpki+Qy/2/'^)Smi,NSai, + l^ai,l/2 + sin {i^ki-Qy /2/'^)^mi,NSai,-lSai,-l/2 

+ COS^(V'fej+Q/^/2/2)(5mi,W<5ai, + 1(5,71, -1/2 + siu^ (V'fci -Q; /2 /2)(5mi ,iV(5ai , - 1 (5^1 4/2 + ^(^ ~ ^l) 

+Sm^,NSa^,-aJau-<y2 [cOs(V'fe_Q^/2/2) sin(V'fe,_Q^/2/2)e-*'3.(fe2-Q„/2)i^-*7^^^_^^_Q^^^^__^^J^^__^^^2 

+ cos(V;fei-Q,/2/2)sin(^,,_Q^/2/2)e^Q^('=^-Q«/2)'o+n4^_^^_Q^j^^__^^_^^__^/^ 
+ cos(V;^^_Q,/2/2)sin(V>^^_Q,/,/2)e-'3i('=-'3:./2)'o-'^'4i,.2-Q;'5ai,+i<5<7i,-i/2 

+ cos(Vi;,_Q, /2/2) sin(V>^^_Q, /2/2)e^«i('=^-'3^/')'°+^^'<5fe2,fei-Q;,^ai,-i<5.i,i/2] , (81) 

which is basically the same as Eq. (64) for a wide well system a.t u = 2N + 2, except for the different quantum 
numbers associated with the Landau levels. 



A. Hcirtree-Fock variational energy 



The full Hartree-Fock energies (noninteracting, Hartree, and Fock energies respectively) can be obtained as follows: 
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and 



E^^ = l(-AsAS + a;.)(ei(0) - 62(0)) + J:(Asas + a;.)(e'i(0) - e'2(0)), 



+2cos{Qyqjl)Qi{q,,)Q2{qn)+'2cos{Q'yqJl)Q'^{qn)Q'2{qn) 

+2cos((0y - Q'y)qnil/2) (ei(g„)e'i(q„) + e2(<z„)e^(q„)) 

+2cos((0^ + Q'y)qJll2) (ei(g„)e^(g„) + Q'M^2{qn))\ , 



E?' = 



-1 



2j^^ ^V^ff^^Niq±)'^cos{qnqyll) &i{qnf + &2{qnf + 2cos{{Q^qy - (5j/9x)^o)Q3(9n)' 

q± n 

+ e'lfeO' + ^'^{qnf + 2cosiiQ',qy - Q'yq^)fM{qnf 

_2 

+ 90~ XI ^JVJV,JViv(9J-) ^ COs(9„gyZo) C0S(9„((3;/ + (3y)Zo) [0l(3n)S2(3n) + B2(5n)0i(9n)] 



+ 



- X] ''^JVJV,JVJv(9l) C0S(7 + 7') ^ e3(g„ - (5x)03(3n + Q'x) COs((feQi, - an5j/)^o)i 



9-L 



(82) 



(83) 



(84) 



where the last term of the exchange energy plays the same role as the last term of Eq. (67) in the wide well system 

d± V = 2N + 2, coupling the nearest and the next nearest pairs of levels to stabilize a coherent phase. The minimum 
of the HF energy is achieved for 7 + 7' = 2rmr {m is an integer) and arbitrary 7 — 7'. 



B. Nonstripe phases 

For the sake of completeness and later discussion, here we first briefly review some results of the CAF phase in 
a double well system in the presence of in-plane magnetic field at even filling factors in our theory. We note that 
this system has also been studied recently in Ref. [36] by numerically solving a single electron HF equation. As has 
been mentioned in Section IIIB 1, in the energy eigenstate basis we use in Eq. (14), the in-plane magnetic field has 
been incorporated in the phase difference between the electrons in the right layer and in the left layer, so that a 
fully incommensurate state [36] has been automatically excluded in our trial wavefunction (see discussion in Section 
IIIBl). The HF energy of such a fully incommensurate state, however, can still be calculated if we let tunneling 
and in-plane magnetic field to be zero in the above HF energy expression and let P± = Q± = Q'j^ = 0, since it is 
well-known that the tunneling energy is effectively absent when the in-plane field exceeds the critical value for the 
commensurate-to-incommensurate phase transition. The HF energy of a fully commensurate state is calculated by 
setting Q± = —Q'^ = in the HF variational energy, and the energy of a partially commensurate/incommensurate 
state is obtained by setting ^ —Q'± or = —Q'j_ ^ 0. The stripe phase functions, tjjk and tl)'^, are just constant 
variational parameters for studying the nonstripe phases. 

Inspecting the last term of Eq. (84)), it is easy to see that the minimum energy should be always at Q± = —Q'^ 
to stabilize the CAF phase even in the absence of in-plane field, i.e. we need only consider Q'^ = —Q± for both 
fully commensurate {Q±_ = 0) and partially commensurate/incommensurate {Q±_ ^ 0) states. (In the terminology 
defined in Section HI A 2 of this paper, these are just the isospin coherent and isospin spiral phases respectively in 
the noninteracting energy eigenstate basis. As stated above, we are using the conventional terminology for most of 
our discussion in this section to avoid confusion.) Therefore, focusing on the nonstripe phase, the HF energy shown 
above can be simplified to be 

Ehf"{Q±) = i^sAs - uj,)Q2{0) + {AsAs + ^.)e'2(0) 

+IeUo, 0) [em' + 02(0)2 ^ Q, (0)2 ^ Q, (0)2] ^ E'AQ±,o) [63(0)^ + Q',{of] 
+^^^(o.o)[ei(o)eJ,(o) + e2(o)e'i(o)] + 2£;^(Qi, 0)93(0)6^(0), (85) 

where the analytic expression for E^^°p{Q±, 0) is shown in Appendix F. From Eq. (85) we find that, as expected, the 
Hartree energy does not contribute to the HF energy of a uniform phase, and only exchange energy (negative value) 
is relevant. 
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By comparing the HF energy of the three different phases: fully incommensurate (A<j/i<j = and P± = Q± = 
0), fully commensurate {Q± = 0), and partially commensurate/incommensurate {Q± ^ 0) phases, our numerical 
calculation shows that the ground state is always the fully commensurate phase (which has the lowest energy), 
i.e. Q± = Q'j^ = as an isospin coherent phase in the noninteracting cigenstate basis. Therefore wc do not 
find the commensurate-incommensurate phase transition (as observed in the odd filling systems) in this even filling 
factor situation. Following the arguments in Section II C (see particularly Eqs. (43)-(45)), we find that this fully 
commensurate state breaks both parity and spin rotational symmetries (but not the translational symmetry). To 
understand that these are two separate broken symmetries one can consider a non-HF state that has uniform ipk and 
V'fc, and (7,7') fluctuating between (—9/2,6/2) and (tt — 0/2, it + 9/2) with fixed 9. Such a state clearly breaks the 
spin symmetry, but not parity (which changes 7 and 7' by tt, see Eq. (43)), and may be described as a spin nematic 
state. 

The equation for the phase boundary between the fully commensurate phase (i.e. commensurate CAF phase) and 
the symmetric phase (both filled levels are orbital symmetric, a = +, but with different spin directions) is 

ojA^sas) = ^(AsAS-S°(0,0))2-£° (0,0)2, (86) 

where cjz is the Zeeman energy. Analogously the equation for the phase boundary between the commensurate CAF 
phase and the fully spin polarized state is 

cofiAsAs) = /b°^((W + A|^ - EUO, 0). (87) 

The calculated phase diagrams of different filling factors and different in-plane magnetic fields are shown in Fig. 13. 

We can also use a similar method to calculate the triplet magnetoplasmon modes by taking 82(0) and 63(0) to 
be small in Eq. (85), because the noninteracting Hamiltonian has been diagonalizcd in the energy cigenstate basis. 
However, since such calculations already exist in the literature both for zero in-plane magnetic field and in the presence 
of an in- plane field [40], we do not show the calculated magnetoplasmon dispersion in this paper. 



C. Exploration of the stripe formation 

When considering the possible stripe formation via the periodic functions of ipk and in the trial wavefunction, 
we need to retain the Hartree energy, which is canceled in the uniform non-stripe phase discussed above. Assuming 
Q± = as in the uniform phase the Hartree-Fock energy in (Eqs. (82)-(84)) for the striped case becomes 

Ehf^O) = i^sAs - u;,)Q2iO) + {AsAs + a;.)e'2(0) 

+1 E ^HiQn) [QM + e2(g„) + e'M + QMf 

+ 2E^(0,g„) [ei(g„)e'2(g„) + 02(g„)e'i((7„) + 2e3(<z„)e^(g„)]} , (88) 

where — 2'Kn/a and a is the stripe period. We define Elj{qn) = {'^T^lQ)~^VNN,NNiln) and Ep{qn) = 
^^^J2q^ ^NN NNi*!-^) ^^^ilnly^o)' neglcct the divergent Hartree energy at qn = 0, which is canceled by the 
background positive charge (providing the overall charge neutrality). Using the perturbation method developed in 
Section HI C, which is equivalent to studying the mode softening of the Goldstone mode inside the CAF phase, we 
obtain the following perturbation energy matrix (see Eq. (56)): 



{AsAS - t^z) cos(V'S) + [eUO) - Ej,{q)] - E°piO) cos(V'o* + Vo*) 
-E"p{q)cos{ro+%*) 

-E°piq) COs(V'o*+^^*) 
{AsAS + ^z) cos(V^,*) + [i?^(0) - E'p{q)] - E-p{0) cos(Vo* + ^^,*) 



where iJjq and ^Jq* are the optimal values obtained from minimizing the total energy inside the commensurate canted 
phase region, and q is the wavevector of the test small stripe as shown in Eq. (50). For simplicity, we have assumed 
that the stripe periods of the two stripe phase functions, tpk and ip'/., are the same, 2n/q. The sign of the eigenvalues 
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of Eq. (89) then determines the existence of stripe formation as discussed in Section IIIC. Our numerical calculation 
shows that both of the eigenvalues of EpJ^f{q) in Eq. (89) are always positive within the canted phase region, showing 
that no stripe phase should occurs with an energy lower than the commensurate CAF phase. We have also studied 
the possibility of first order phase transition to a stripe phase, which cannot be included in the perturbation theory, 
by directly comparing various ground state energies. We still find that no stripe phase can be stabilized energetically. 
Therefore, unlike the rich phase diagram shown in the odd filling factor double well system at = AN + 1, the system 
at even filling v = 4iV + 2 in the presence of in-plane field has neither the commensurate-to-incommensurate phase 
transition in the large B|| region (see discussion below Eq. (85)) nor the uniform-to-stripe phase transition in the 
large layer separation region. The former result is also consistent with the recent result obtained in an effective field 
theory [39] . We discuss this result further in Section VIII. 

VIII. DISCUSSION 

In this section, we discuss and compare the quantum phases obtained in different systems studied in this paper (as 
well as those already existing in the literature). 

A. Directions of isospin stripe and isospin spiral orders 

In Table II we summarize the isospin spiral and isospin stripe directions obtained in the three different cases, Wl' , 
W2' , and Dl, where at least a spiral or a stripe order exists in our HF calculation. For the convenience of discussion, 
we will define the isospin components to be in the layer index basis for the double well systems ai v = AN + 1, i.e. 
the isospin spiral/coherent (stripe) phases in Dl case correspond to the conventional commensurate/incommensurate 
(stripe) states respectively. We do not include the commensurate CAF phase in D2 system, because it also has a 
spiral order following the in-plane field in the layer index basis and therefore is similar to the commensurate state 
(isospin spiral phase) of Dl system. 

Note that although the definition of isospin components are different for these systems (see Table I), all the isospin 
spiral orders select a wavcvector, Q± (but ±Q±_ are degenc;ratc; for the wide well system), to be perpendicular to the 
in-plane magnetic field. The mechanisms for the spiral winding in wide wells and in double wells are very different. 
In a wide well system with a strong in-plane magnetic field, the electron wavefunction is distorted by the anisotropic 
field and therefore the electron-hole binding energy is the strongest if the isospin is winding along the direction 
perpendicular to the in-plane field [17]. More precisely, this is true only in the presence of a very strong perpendicular 
magnetic field or in wells which are not too wide, so that the effective cyclotron resonance energy, ui±, is comparable 
to or larger than the confinement energy, luq. The direction of the spiral winding can change to the in-plane field 
direction if the perpendicular field is so weak that loj_ <^ Uq- (We will show this result explicitly in below). This 
effect of spiral winding locking is less important in the weak region, leading to the isotropic; phases observed in 
the intersubband level crossing region (VFl). In a double well system with zero well width, however, the in-plane 
magnetic field does not affect the single electron wavefunction, but only affects the tunneling amplitude through the 
Aharonov-Bohm phase factor, which selects a specific spiral wavevector, Q± = {0,Px) = (0, ec?i3||/c), perpendicular 
to the in-plane magnetic field. When B|| is weak, the isospin spiral order follows the wavevector of the tunneling 
amplitude, since it minimizes the tunneling energy and does not cost much in the Coulomb exchange energy. When 
_Bj| exceeds a certain critical value, the energetic cost of winding from the point of view of exchange energy becomes 
prohibitively high, and the system goes into the incommensurate phase with no isospin winding. Therefore the origin 
of isospin spiral order in a wide well system and a double well system is very different, although both are perpendicular 
to the in-plane magnetic field. 

Now we discuss the mechanism, which determines the directions of the isospin stripe orders (if it exists) in these 
systems. From Table II wc note that the stripe modulation of Iz (i.e. the normal vector, n) of the stripe is along x 
axis, parallel to the field for a wide well system at both even and odd filling factors (for even filling system, W2', we 
cannot stabilize a many-body phase within the HF approximation, but speculate, based on very small calculated HF 
energy diffcrcnc;cs, that the resistance anisotropy observed in the experiment [15] could result from a skyrmion stripe 
phase near the level crossing region, which could perhaps be stabilized by going beyond our approximation scheme 
(see Section V). But the stripe normal vector is along y axis, perpendicular to the By field, for a double well system 
aX V = AN + 1 (no stripes are found in a double well system sX v = AN + 2). The different stripe directions in the 
two systems could be understood as a result of competition between two effects: one is the anisotropy energy (i.e. the 
total energy difference between a stripe perpendicular to the in-plane field and a stripe parallel to the in-plane field, 
see [11]) induced by the in-plane magnetic field via the anisotropic distortion of electron wavefunction, and the other 
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one is the exchange interaction between the spiral order (oscillation of {I^) and {ly)) and the stripe order (oscillation 
of (/z))- This effect reflects on the term, S|'^^((?n, ; Q^;, Qy) of Eq. (60) for a wide well at v — 2N + 1 and on the 
last term of Eq. (75) for a double well at v = AN + f (the similar terms in the even filling systems can be also found 
in Eq. (67) and Eq. (84) for wide well and double well systems respective: proportional to ^^{qnY and Q'^{qn)^)- 
Such interaction between the spiral order and the stripe order prefers to keep the wavevectors of these two oscillations 
parallel with each other in order to optimize the exchange energy. On the other hand, the anisotropy energy induced 
by the anisotropic electron wavefunction in the finite width well prefers to form a stripe aligned perpendicular to the 
in-plane field direction, i.e. its normal direction is parallel to x axis [58] . 

If the spiral-stripe exchange energy dominates the anisotropic energy (like in the double well system, where the 
zero- well- width electron wavefunction is isotropic and hence the anisotropy energy is zero) , the stripe order prefers to 
select a normal wavevector parallel to the wavevector of the spiral order, resulting in a spiral stripe phase. The stripe 
direction in this case is governed by the spiral order, which is perpendicular to the in-plane field in the commensurate 
phase as we find in the numerical results shown in Section VI and in Fig. 12(a) (phase (IV) and (V)). When the 
in-plane field is so strong that the double well system undergoes a first order phase transition to an incommensurate 
phase, the spiral order disappears and therefore the stripe direction is not locked by the in-plane field, being in an 
arbitrary direction in the two-dimensional plane (phase (III) in Fig. 12) [24]. On the other hand, when we consider a 
single finite width well in the presence of in-plane field, the anisotropy energy is finite and competes with the spiral- 
stripe exchange energy. If the anisotropy energy dominates, the spiral-stripe exchange energy can not force the stripe 
and spiral orders to be parallel with each other (and hence perpendicular to the in-plane field). Therefore the density 
modulation of the stripe phase prefers to stay in the direction of in- plane field [58] , having a direction perpendicular 
to the spiral order, resulting in a skyrmion stripe phase. Therefore the mechanism we discussed so far can explain all 
of the spiral and stripe directions obtained individually in the previous sections within the HF approximation. 

As a final remark, we note that the HF energy difference between the spiral stripe phase and the skyrmion stripe 
phase are very small compared to the other energy scales in our numerical calculation, both in wide well systems 
and in double well systems. Therefore we expect that various effects (e.g. disorder and impurity scattering) ignored 
in our calculations could have strong influence in determining the eventual ground state of the system, providing an 
experimental stripe direction different from that obtained in our idea HF theory. The actual finite width effects of 
a double quantum well system in a realistic experiment could also lead to the stabilization of the skyrmion stripe 
(rather than a spiral stripe) phase, especially when the spiral-stripe coupling is weakened by the finite temperature 
and/or finite disorder effects. 



B. Role of spin degree of freedom in level crossing 



From the results presented in Sections IV and V for wide well systems, we find that systems of odd filling factors 
iW\ and W\') have more interesting coherent phases with exotic quantum order than the systems of even filling 
factor {W2 and W2'). As summarized in Table I, the HF trial wavefunction, Eq. (36), stabilizes an isospin coherent 
phase in Wl , and an isospin skyrmion stripe phase in WV . For systems at even filling factors, strictly speaking, no 
novel quantum phases are obtained in the wide well systems within our HF approximation (although sometimes the 
exotic phases arc very quantum closeby in energy). Since the orbital wavefunctions of the two crossing levels in the 
even filling systems are the same as those in the odd filling systems, it is natural to attribute the important difference 
between the two systems to the additional spin degree of freedom of the two crossing (or degenerate) levels, that is 
present in the Wl/Wl' cases. More precisely, we can compare the formulae of the HF variational energy of an odd 
filling system shown in Section IV with the formulae of HF energy of even filling system in Section V, where we can 
simply take ■0^ = in Eq. (49) and Eq. (64), and consider the trial wavefunction constructed by the two crossing 
levels only (instead of the four degenerate levels), similar to Eq. (36). In such situations, the real spin quantum 
number of the two isospin states are of opposite sign, different from the level crossing in the odd filling system. The 
resulting HF variational energy (denoted by underlines) then becomes (compared to Eqs. (58)-(60)) 



(90) 
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+^ E "^nM.n. cos((Q,9^ - ^ cos(9„g^/2^e3(g„)^ (92) 

where we have neglected those uniform terms (hnearly proportional to 0^(0) for i = 1, 2, 3) for simplicity. Comparing 
Eqs. (90)-(92) with the HF energy of an odd filling system in Eqs. (58)-(60), we find that the odd filling factor 
systems have one additional term in both Hartree and Fock energies (i.e. the E]^q and Epf^ terms) after neglecting 
the imiform terms and those singular terms proportional to (5q,^,o shown in Eqs. (58)-(60)). These two additional 
terms result from the direct and the exchange contractions of the correlation, (ct c- ct ,c- ,), which is absent 

^ ' \ ni,(7 n2,cr 712, cr' ni,cr /' 

at even filling factors, because the two relevant levels, ni and n2, have opposite spin direction. Since such additional 

exchange energy is larger than the additional Hartree energy, a coherent phase near the level crossing point can be 
stabilized more easily in an odd filling system than in an even filling system. This mechanism explains the results we 
obtain in the wide well system within the HF approximation, although it does not exclude the possibility of having a 
many-body phase in an even filling system in more sophisticated theories. 

Such effects of spin degree of freedom can also be observed from the results of the double well systems {Dl and D2). 
When considering only the two crossing levels in 1/ = 4N + 2 case, we will not obtain any many-body state for the 
same reason as mentioned above. The coherent phase (i.e. commensurate canted antiferromagnetic phase) obtained 
in the literature [33,34,36] and in this paper strictly relies on the incorporation of all the four degenerate levels (see 
Sections HIB and VH). On the other hand, at least four different many-body phases (see discussion in Section VI and 
Fig. 12) are obtained as ground states a.t v = AN + 1 within the HF approximation. In the odd filling case, electrons 
in the top level equally occupy the two layers, and therefore when the layer separation is larger than the order of Iq, 
interlayer coherence becomes relatively weaker and a stripe formation may occiir in order to optimize the intralayer 
exchange energy in each layer, refiecting the charge density wave instability studied in the high half-odd-integer single 
layer systems [4,5] . Such a simple mechanism, however, does not seen to apply to u = AN + 2 case, because every 
flux quantum is occupied by an electron in each layer — any nonuniform density modulation, as in a stripe phase, 
will have to pay a large direct energy for the double occupancy of electrons of opposite spins in each layer in each 
flux quantum. Since the commensurate canted phase can successfully lower the exchange energy by canting both spin 
and isospin degrees of freedom, a stripe phase becomes unlikely in the even fllling factor y = AN + 2 bilayer case (at 
least in the absence of any external bias voltage). When an external gate voltage is applied in the growth direction 
of the two layers, the two layers become partially flUed, and there might be some interesting stripe phase in this 
situation, at least for the following extreme case: when the gate voltage between the two layers is so large that the 
fllling factor of the left layer is near 2 A'' +1/2 and of the right layer is near 2 A'' + 3/2, the two layers can form stripe 
phases individually similar to those well-known single layer systems at high half-odd- integer filling factors {N > 1) 
[1] (especially if the layer separation is so large that the tunneling energy and interlayer interaction are negligible). 



C. Comparison with the exchange induced spin density wave in a very wide well 

It is very interesting and instructive to compare our results of stripe phases with the spin density instabilities 
proposed first by Brey and Halperin (BH) [44] in a wide parabolic well system and by Das Sarma and Tamborenea 
[59] in a zero-field double quantum well system. In Ref. [44], it is found that when an in-plane magnetic field is 
applied (without any perpendicular magnetic field) to a very wide parabolic well (width is 4000 A), Hartree energy 
between electrons inside the well modifies the electron density profile to be almost a uniform slab, making the parabolic 
well a good approximation of the 3D jellium model where electrons move in a constant positive charge background. 
Calculation of spin density correlation function shows that a divergent singularity occurs at a certain wavevector in 
X direction (parallel to the in-plane field), Q^^ = kF,i + fc_F,x, where kF,a is the Fermi wavevector of spin a electrons. 
They proposed that such spin density wave instability is similar to those originally proposed by Celli and Mermin 
[60] in three-dimensional systems and may cause a resistance anisotropy in the x — y plane of the well. This result 
[44] seems to contradict the spiral direction suggested by the magneto-roton minimum of the magnetoplasmon mode 
[17] as well as our HF calculation in this paper for a wide parabolic well system at the even filling factor, where the 
softening of the spin triplet magnetoplasmon mode (i.e. the divergence of a spin density correlation function) is in y 
direction, perpendicular to the in-plane magnetic field. Besides, the wavevector of the BH spin density wave mode 
softening, Q^^ , depends mainly on the electron density and is almost independent of the magnetic field strength, 
while the wavevector of the mode softening in our calculation crucially depends on the in-plane magnetic field strength 
and is almost independent of the electron total density [17]. Here we analyze the superficial contradiction between 
our results in this paper and the BH results, providing a deeper understanding of these two calculations. We find that 
both results are correct, and their different characteristic spin density wave vectors arise entirely from considerations 
of different limits of system parameters in the two situations. 
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First wc should clarify the energy scales of the wide well system we are going to discuss. For the quantum Hall 
situation of our interest in this paper, the perpendicular magnetic field is always finite {B±^ ^ 3 Tcsla, i.e. w_l ~ 5 meV) 
to ensure the existence of a discrete orbital Landau level spectrum. The well width is chosen to have a comparable 
{u>o ~ u>±) bare confinement energy ~ 7 meV. The (iso)spin density wave instability (for both odd and even filling 
factor) in this system occurs in a very high in-plane magnetic field region ~ 10 — 30 T, or ~ 17—50 meV). 
Therefore we can obtain the following approximate order of energy scales for the system of our interest in the current 
paper (using the even filling factor system as an example) 



u± u)o <^ and Ni ^(1)) 



(93) 



where is the Landau level index of the highest filled level of spin a, and the total filling factor ly = + + 2 
is of order of unity, A small parameter, e = Maa;(a;_L, wo)/w|| *C 1 can therefore be defined. In this limit, we have 
oji and e^t^^ji for the two noninteracting Landau energy separation (see Eq. (2)). Besides, the two effective 

magnetic lengths, li^2 = -\/l/m*a;i,2, defined in Eq. (5) (or see Ref. [17]) become Zoe^^^ and /qc"^^^ respectively. On 
the other hand, for the system discussed by BH [44], no perpendicular magnetic field is applied and the predicted spin 
density wave instability occurs at an intermediate strength of 1 Tcsla, i.e. ~ 1.7 meV) for a very wide well 

(for width about 4000 A, the parabolic confinement energy uq < 0.05 meV). Therefore we can obtain the following 
approximate order of energy scales for the BH system: 



a;_L <C Wo < <^||, and ATt NtotA, N\ iVj 



(94) 



where A^tot.Kj.) is the total electron number of spin up(down) (in the direction parallel to B||), and hence the Landau 

level degeneracy, A^, is only of the order of unity in this weak B± limit. Therefore there are two small parameters 
we can construct for the BH system: e = Max(u)±,uJo)/u)\\ = uJo/^^\\ <C 1 as defined earlier and a = (wx/wq)^ <C 1, 
which is of order unity in the system of our interest in this paper. In this limit (Eq. (94)), 0)1^2 are still the same 
as above (i.e. close to a;|| and e^cjn respectively), while li lae^^^a and I2 l^e^^^'^a. Note that although both 
e and a arc small numbers for the BH system, an additional constraint, a <C e, exists also, to ensure the fact that 
no perpendicular magnetic field is applied. We will show below that this additional constraint (i.e. a <C e) and the 
resulting small Landau level degeneracy {N^ ~ C'(l)) are the key points needed to understand BH's result of SDW 
instability from the perspective of our theory. 

For the convenience of comparison, we still consider the two level degeneracy on the top of the filled Landau levels 
and use the "core state" approximation for a qualitative discussion, although it is certainly a bad approximation in 
the limit of Eq. (94), where the Landau level energy separation goes to zero. To simplify our analysis, we calculate 
i?]^2((5x) shown in Eq. (69) only, which corresponds to the electron- hole binding energy via ladder diagrams and 
gives rise to the roton minimum of the magnctoplasmon dispersion [17], and use a zero-range contact electron-electron 
interaction. We can obtain the following approximate E^'^{Q±) by using the analytical expression of A{q) in Eq. 
(Bl) and taking the e — > limit: 
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where we have let n\ = (0, A^j^) and 712 = (0, A'^ -|- 1) for the orbital level index of the top filled and the lowest empty 
levels respectively. 

For the quantum Hall system of our interest in this paper, a is of the order of unity (see Eq. (93)), and then the 



above equation can be simplified further (denoted by E^^'^{Q±)) by keeping only the leading term in e: 
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where we have used oc x" ior x ^ 1, = Ni = N = (i^ — 2)/2 for even fiUing system, and I2/I0 ^ ct/y/^- It is 

very easy to see from Eq. (96) that E^^'^{Q±) is nonzero only at Qa; = and a finite Qy. The extreme value of Qy 
to optimize the exchange energy is given by the length scale generated by the competition between the exponential 
function and the power-law function, proportional to \/Nl2^. This is consistent with our numerical results obtained 
either from the HF variational calculation in this paper or from the collective mode calculation in our earlier work 
[17]. 

On the other hand, if we take a ^ first in Eq. (95) as appropriate for the situation considered by BH, we obtain 
the electron-hole binding energy (denoted by E^^'^{Q^)) to be 



Ex {Q±) ~ -o-Z^exp -— exp — ^ L^^+i L^^ cos{{q^Qy - qyQ^)lo) 
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^B- 5q^,oC / ^ cosi./m'^qyll/h) cosi^/Wlqyll/h) cos{qyQJl), (97) 

Jh/il %'o 

where we have used the asymptotic formulae, (a;) oc e^/'^x~^/^ cos{2^/nx) for large n [61], and have let the constant 
contribution from integration of small qy region to be denoted by B; C is a constant factor. It is easy to see 
that the max;imum electron-hole (exciton) binding energy occurs at Qy = and ~ ±(fci?.-f -I- kp^i), where the 

Fermi wavevector is determined by the Fermi energy, kF,a = \/2m*{EF,a — Ek=o) ~ ^2m*{E^^ — E'^^^ ^) = 

V2m*N„uj2 = \/2iV^/^\ for spin (J subband {Ek = k"^ /2m* is the usual free particle zero-field energy dispersion with 
effective mass m*). Therefore Eq. (97) is consistent with the result obtained by BH, showing a SDW instability along 
the in-plane field direction (x) with a wavevector, Qx = ±{kp_if + kF,i)- Note that the maximum E^'^''^{Qx,0) does 
not occur at = i(^F,t ~ kF,i) because the decaying function, q~^ in Eq. (97) has higher contribution from the 
small qy region, and hence a larger value of \Qx\ (i-e. \kF,i + fc-F,il) should a give stronger electron-hole binding energy 
than a smaller value of \Qx\ (i.e. jfe^,! — kF,i\)- This effect is enhanced for the long-ranged Coulomb interaction, 
which has a faster decaying function in qy. Therefore our simple analysis above shows that the apparent inconsistency 
between of the SDW spiral winding direction studied in our work [17] and in the BH work [44] is just due to the 
different limits of interest in the two works: A quantum Hall situation with finite B± in our case and the "zero-field" 
B±=0 situation in Ref. [44]. 

We now summarize the above discussion on the spiral (spin-density- wave) order of a wide well system in the presence 
of in-planc magnetic field (along x axis) : when a strong perpendicular magnetic field is applied so that only few Landau 
levels are occupied and the in-plane field is tuned to be close to a level crossing (or a level near degeneracy) point, 
the (iso)spin spiral winding prefers a specific wavevector in a direction (y) perpendicular to the in-plane field, while 
the (iso)spin (7^) oscillates in x direction; when the perpendicular magnetic field is reduced to a very small value, the 
electron wavefunction near Fermi energy becomes a plane wave and strongly modifies the interaction matrix elements 
through the form function — the spin-flip exchange energy then stabilizes the spin density wave at a wavevector 
Q* ~ ±{kK,-[ + kF,i)- The spin (or isospin) density wave state actually results from two different mechanisms in the 
two different limits of the same system. These results are independent of the gauge choice in the theory. 



D. Effective Hamiltonian of wide well at u = 2N + 1 ciround intersubband level crossing point {Wl Ccise) 

In the previous sections, we have been using the concept of "isospin" to discuss the ground state wavefunction 
properties in both the wide well systems and the double well systems. Without losing generality, one can also apply 
the isospin concept to effectively represent the Hamiltonians of these systems so that some aspects of this physical 
properties (e.g. collective excitations or finite temperature effects) may be studied in a more comprehensive way. 
Here we will give the effective Hamiltonian of a wide well system a.t v = 2N + 1 near the intersubband level crossing 
point in the isospin representation after projecting the whole system into the two crossing levels we consider in this 
paper. According to our HF calculation results presented in Section IV, only the uniform coherent phase is stabilized 
in this intersubband level crossing region, and therefore we will not consider the stripe order (modulation of isospin 
z component) here for simplicity. 

For convenience of comparison, we redefine the isospin component to be (c.f. Eqs. (38) and (39)): 
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where a = x,y, z and crf_^j^ is the Pauli matrix clement. Therefore the low energy effective Hamiltonian for a wide 
well system at v = 2N + 1 can be obtained from Eq. (9) by using Eq. (98) and considering the long wavelength limit; 

HY/f = J dr^ |(/x + 5.o)^r +5.. (^r)' + +5xCOs(27) +5a/3(7)(Va7)(V;37)] {iY')'] (99) 

where we have used — cos 7 and — ^^^1 to address the phase dependence in the isospin x — y 
components. Va is the spatial derivative of phase 7 in a = a;, y direction. The coefficients shown above are respectively 
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where we have used the fact that "l^^yg^ f^^i (0) and f|-4j.(0) {^^a is the derivative in the a = x{y) component of 

wavevector, q±) are pure real, and dxdyV^^^^{0) is pure imaginary giving an additional phase 7r/2. The effective 
interaction matrix element, Vj^j^ 7^7^ ((fj_), can be obtained within the Hartree-Fock theory: 

VZ.JsiAQ^) ^ ^VZ.j.iM^) - E ViYiMiP^) e'^^^^^-^'^y^^l. (106) 

From the effective Hamiltonian in Eq. (99), we can see that the effective "Zeeman" energy, h + Qzq, is the same as 
the single electron energy gap between the two crossing levels (including the HF correction). Therefore when such 
"Zeeman" energy is dominant in Eq. (99) (i.e. the two levels are far from crossing region), the isospin direction must 
be polarized at = ±1/2 and hence the ground state has no transverse component, i.e. lY^ = 0. These are the 
conventional integer QH states with electrons of the top filled Landau level being at level n = (1,0) (for isospin up) 
and at n = (0, TV) (for isospin down) respectively. In the intermediate region between the two states (sec also Fig. 
7(b)), where \h + gzo] is small compared with the other energy scales, the transverse isospin component, J-Y^ rnay 
become finite to minimize the energy so that the isospin polarization becomes canted. The prefactor, cos(27), of the 
g-^^ term in Eq. (99) selects 7 = 0, tt to be the two degenerate points of the lowest energy, showing that only lY^ 
is finite and there is no isospin in y direction. This two-point degeneracy results from the fact that only the parity 
symmetry is broken in the intermediate region. The coefficient of the gradient term, Qapil), is positive in such an 
intersubband level crossing case so that the optimal value of phase 7 is always a constant in space, confirming the 
fact that no spiral order exists in the coherent phase. 

For the convenience of comparison, here we also show the effective Hamiltonian of the double well system at 
y = 4A^ + 1 using the same definition of isospins in the layer index basis {B\\ is along x axis): 

n^lf = J dr^ |-A5Ascos(7- Px2/)lf ^ + [V{;,^^(0) - V^^^,^^(0)] (xf 1)' 

+ [2Vf^,4^(0) + dlV^^^^^mV^^f] (Jf , (107) 

where Vj^/^ i^i^iqi.) is similarly defined as Eq. (106) now for the double well system. Without losing generality, the 
effective Hamiltonian of the wide well system in Eq. (99) is similar to that of the double well system in Eq. (107) [26], 
except that in the former case the in-plane field is included implicitly in the interaction matrix element, V/^^ i-^iS^^') 
and the isospin stiffnesses in x and y directions depend on the isospin direction. The effective "tunneling" amplitude, 
g-^ , in the wide well system originates from the Coulomb exchange energy and as such is not an independent physical 
quantity (in contrast to) the tunneling energy in the double well system, which is an independent physical parameter. 
When considering the intersubband level crossing in the small in- plane magnetic field range (i.e. Wl case), gapil) 
is always positive as mentioned above, so that the phase is fixed to be or tt similar to the commensurate phase 
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in the double well system at v = AN + 1, where only 7 = is chosen by the tunneling energy, A^^i^. However, 
if we apply Eq. (99) to the intrasubband level near degeneracy case in the strong in-plane field region (i.e. Wl' 
case), 50/3(7) may become negative, so that the total energy can be minimized at some finite isospin winding in 
the isospin x — y plane (say, j = PxU becomes a function of spatial coordinate). This result leads to the finite 
spiral order of the isospin skyrmion stripe phase as discussed in Section VI in the large in-plane field region. In 
such a case, the effective "tunneling" amplitude, gx, becomes zero since the space average of 003(27) = cos{2Pxy) 
vanishes. Phenomenologically this situation is similar to the incommensurate state of the double well system, where 
the tunneling energy is dominated by the Coulomb exchange energy so that the phase is not fixed by the in-plane 
field. But in such intrasubband level near degeneracy case, we will also have to consider the stripe order (gradient 
term oi T^^) to obtain the a full effective Hamiltonian in describing the skyrmion stripe phase obtained in the HF 
calculation. Finally we note that in general, the interlevel(layer) scattering amplitude (proportional to V^^^^^{q±)) 
is finite in the wide well system, while it is zero in the double well system due to the layer separation. This gives the 
finite values of gx and the second term of g^pi'y)-, according to Eqs. (104)-(105), and leads to the major difference 
between a wide well system and a double well system. 

From the experimental point of view, however, a smooth crossover from a wide well system to a double well system 
can be observed in the same semiconductor system by tuning some parameters, e.g. electron density [62] due to the 
Coulomb screening effects on the confinement potential. It is therefore very interesting to investigate such crossover 
from monolayer to bilayer systems from a more fundamental theory they are usually assumed to be two different 
systems when studied in the literature (the inclusion of a modified form function to take into account the finite width 
effect in the double well system does not help in this comparison). The effective Hamiltonians shown above could be 
a good starting point for this problem, but we will not discuss it further here since it must incorporate one-electron 
self-consistent potential [45,59] in the confinement potential which is beyond the scope of our approximation scheme. 
We must emphasize that it is useful to remember that the bilayer and the wide well systems may be continuously 
tuned into each other by changing the system carrier density. 

E. Instability of the stripe phases and the domain wall formation 

In a single well system at partial filling factors, another interesting kind of charge density wave ground state, the 
so-called "bubble phase", may occur when the filling factor of the top level is away from the exact half odd integer 
values [1]. It is believed that when the filling factor of the top Landau level is below about 0.4, electrons in the top 
level may lower their potential energy by accumulating to become bubbles, rather than stripes, in the two dimensional 
plane [63]. The equivalent "bubbles of holes" may occur when the filling factor of the top level is between 0.6 and 
1. This phase can be understood as the edge state instability of the stripe phase due to the backward scattering 
between the two nearest edges, leading to the breaking up of the stripes into segments. These segments of stripes then 
rearrange their shapes and positions to form bubbles in a lattice to lower their total energy. (The strong field Wigner 
crystal phase can be viewed as an extreme limit of the bubble phase — having only one electron in each "bubble" .) 
A possible experimental evidence for such an interesting phase is the reentrant integer quantum Hall effect observed 
recently [1,64]. Similar nonstripe density modulation may also occur in the integral quantum Hall systems discussed 
in this paper. Near the level degeneracy region, the edge state instability may bo strong enough to break the stripe 
phase and form several "bubble-like" domain walls inside which the isospin (J^) is different from that outside the 
domain wall. Some spiral winding structure may appear at the surface of the domain wall to reduce the Hartree 
energy (see Fig. 14). Such domain wall structures could be stabilized by the presence of surface disorder as suggested 
in Ref. [53], and may also lead to the observed resistance anisotropy [15,16]. However, this domain wall phase cannot 
be included in the trial wavefunctions we propose in this paper [65], and we speculate that they may be important in 
understanding some other experimental results, such as the hysteretic behavior in the magnetoresistance experiments 
[51,52]. 

F. The v = 1 bilayer system 

There has been a great deal of recent experimental and theoretical interest in the (weak-tunneling) bilayer double 
well system at u = 1, i.e. N = case, which has essentially been left out of our consideration in this paper where 
we have mostly restricted om-selves to the case > 1, with = 4N + 1 (Dl) or 4N + 2 (-D2), level crossing bilayer 
situation. (The 1^ = 1 bilayer case can be thought of as an isospin level crossing situation for = in our notation.) 
Here we provide some brief comments on the v = 1 bilayer situation in the context of our theoretical results presented 
in this paper. The specific issue we addresses is the nature of the phase transition in 1/ = 1 bilayer system as the 
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layer separation (rf) is increased in the zero (or weak) tunneling situation. It is commonly believed, not however 
based on any really compelling evidence, that the v = 1 bilayer system undergoes a first order transition from an 
incompressible interlayer-coherent quantum hall state (presumably a Halperin (1,1,1) state [23]) to a compressible 
state (presumably two decoupled f = 1/2 single-layer states) as d increases above a critical value (which depends 
on the tunneling strength). We want to suggest here another distinct possibility based on our results presented in 
this paper. It is in principle, possible for the transition to be a second order quantum phase transition (rather than 
a direct first order transition) to a many-body state with exotic quantum order, such as an isospin stripe phase 
discussed in this paper. Such an exotic phase will still be incompressible (for low value of disorder), but perhaps with 
a much smaller value of the incompressibility gap. (The system will make a transition to the compressible state of two 
decoupled f = 1/2 layers at some still larger nonuniversal value of d.) We believe that there is already some evidence 
supporting our suggested "double-transition" scenario (a second order transition to a weakly incompressible stripe 
phase at d = dd, followed by a first order transition to the compressible decoupled // = 1/2 layers at d = dc2 > del)- 
In particular, the early calculation of Fertig [56] finding a finite wavevcctor magnctoplasmon mode softening in bilayer 
TDHF theory indicates the obvious possibility of an isospin stripe formation at the characteristic wavevector of mode 
softening. The second order quantum phase transition associated with this mode softening leads to an incompressible 
state (at low disorder) which has a finite quasiparticle gap (perhaps reduced from that in the uniform (1,1,1) phase). 
At some higher value of disorder the stripes would eventually be pinned independently in each layer, leading to an 
incompressible-to-compressible phase transition as a function of increasing disorder. (This is also consistent with 
experiment where all the interesting results are typically obtained in extremely high mobility bilayer samples.) The 
level crossing HP technique used in our current work is unfortunately unsuitable to investigate the ly = 1 bilayer 
case, which has to be studied by the direct numerical diagonalization technique. It is therefore encouraging that 
there is recent numerical evidence [66] in apparent support of the scenario as proposed here. The recent experimental 
results of inelastic light scattering [67], showing that the magneto- roton excitation is softened at a layer separation 
very close to the incompressible-to-compressible phase boundary, is also consistent with our scenario that there is an 
intermediate second order phase transition to a stripe phase before the system undergoes the first order transition 
to a compressible state. The magnetoresistancc anisotropy measurement around the critical layer separation of the 
observed magneto-roton excitation softening (which is also very close to the incompressible-to-compressible phase 
transition boundary) should be a good test of such charge density wave order in the v = 1 incompressible quantum 
Hall system. 



IX. SUMMARY 



In this paper, we use the mean-field HP approximation to systematically study the possible spontaneous breaking 
of parity, spin, and translational symmetries in different integer quantum Hall systems in the presence of an in-plane 
magnetic field, concentrating on the level crossing situation in wide parabolic well and bilayer double well systems 
for both even and odd filling factors. We propose a general class of variational wavcfunctions to include the isospin 
spiral and isospin stripe orders simultaneously, and discuss the symmetry breaking as well as the exotic quantum 
order properties of various many-body phases generated by our wavcfunctions. Comparing the HP energies of these 
many-body phases, we find several of them can be stabilized near the level crossing or level near degeneracy regions 
of different systems, breaking certain system symmetries as listed in the last two columns of Table I: (i) for a wide 
well system at v = 2N -|- 1, we find an isospin coherent phase in the intersubband level crossing region (small S||), 
breaking the parity symmetry only, and an isospin skyrmion stripe phase in the intrasubband level near degeneracy 
region (large B\\), breaking the parity and the translational symmetries in both x and y directions simultaneously 
and having finite topological isospin density, which leads to appreciable charge oscillations in the direction parallel 
to the in-plane field; (ii) for a wide well system at v = 2N -|- 2, we find direct first order phase transitions between 
simple (un)polarized QH states in both intersubband and intrasubband level crossing regions, but we suggest that 
the resistance anisotropy observed recently by Pan et. al. [15] may possibly be explained by a skyrmion stripe phase 
by going beyond the HP approximation; (iii) for a double well system at u = AN + 1, we stabilize the coherent, 
spiral, coherent stripe, and spiral stripe phases in different parameter regions of the phase diagram (see Fig. 12), and 
critically discuss the broken symmetries and the exotic quantum order of these many-body phases in Section VI D; 
(iv) for a double well system at AN + 2, only a coherent phase (= commensurate CAP phase in the literature) is 
stabilized, breaking parity and spin rotational symmetries simultaneously. We also compare our HP results for these 
different systems in details, and discuss the influence of wavefunction anisotropy, spiral-stripe coupling, spin degree 
of freedom and finite well width in Section VIII, manifesting the nontrivial effects of Coulomb interactions in the 
multi-component quantum Hall systems. 
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APPENDIX A: THE VARIATIONAL WAVEFUNCTION FOR STRIPES PARALLEL TO THE IN-PLANE 

FIELD 

In this section we derive a variational trial wavefunction, which describes a stripe phase parallel to the in-plane 
field in a wide well system. As mentioned in Section III, we choose the Landau gauge, = (— _Bxy, — _B||Z, 0), 

in which particle momentum is conserved along x axis. We will use the notation "O" to denote quantities, O, in this 
gauge in order to avoid confusion. The noninteracting Hamiltonian, Eq. (6), becomes 



1 / eB^yV 1 / 



eB\\z 



Pz 



1 



c J 2m* 2 



(Al) 



Setting = to be a good quantum number, we can obtain the single electron wavefunction similar to the form of 
Eq. (4): 



<i>^{y-t',k,z), 



(A2) 



where ^fi{y,z) satisfies the following wave equation: 



1 » 2 2 , ^ f , 1 * 2 2 c 

-m u , y H py — H — + -m ojnZ — 

2 2m* c 2m* 2 ° 



^^{y,z) = E°^f^Wiy,z). 



(A3) 



If we define an auxiliary function, ^^'{pIq, z) to be 



dy^^{y,z)e-^py, 



(A4) 



then it is easy to show that '^'^'(pIq,z) satisfies exactly the same equation as '^^{x,z) in Eq. (1) by redefining 
X = X + Mq. Therefore we could write down the complete solution of Eq. (A3) to be 



{y, z) = ^j dp {pll cose-z sin 6) ■ {pl^ sinO + z cos 0) , 



(A5) 



where rpn\x) has been defined in Eq. (5). Note that the energy quantum number, n, is the same as before because 
the energy eigenvalue is independent of the gauge we choose. 

Similar to our analysis in Section III A, we can construct the following trial wavefunction for the possible many-body 
ground state near the degeneracy point of the two crossing levels (or of the two nearly degenerate levels): 



g^ 

C2,fe 



^-ikQyil/2-i^/2 cos(V/c/2) e''='3'''o/2+'^/2 sin(V'fe/2) 

_g-ikQyll/2-ij/2 sin(V;fe/2) e''='3«'o/2+i7/2 cos(Vife/2) 



5^ - 

|,fe-Qx/2 



Then the isospin components can be obtained easily in the new ground state wavefunction: 

{Uq±)) = In^S,^,o te(g-l,0)e'«»'3^'o/2e,(^,)- A|r.(5-l,0)e-*«^«^'o/2e,(^^) 



(A6) 



(A7) 



and 
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(J+(g-l)) = 7V^^,^,Q^A|^^(g-l, 0)63(9, - Qy) e'^. (A8) 

Therefore the stripe phase constructed in this method is along x axis and the spiral direction is determined by Qj^ . 

Finally it is instrTictivc to point out that the form function obtained from the noninteracting electron wavcfunction 
in Eq. (7) is invariant under such gauge transformation. Defining the form function in the new gauge similar to Eq. 
(7), we have 

= jdxj dze-"'''^e-"'^^^l{x-llqy/2,z)<^l{x + llqy/2,z) = AW^^{q), (A9) 

where we have changed the integration variables: p = qy/2 + x/Iq and p' = qy/2 — x/Iq in the last equation. Therefore 
the interaction matrix element, V^ft^^fi^fi^iq) does not change in the new gauge, and the HF variational energies 
for the wide well system in the new gauge, A[^^{f), can be obtained by doing the following simple transformation 
in Eqs. (59), (60), (66), and (67): V^.^.^.^{q^,0) V^^,^ ,^^^(0, g„), V^.^.^.^{q^,Qy) V^^mM^^-)' 
Qsiqn-Qx) — > Q^iqn-Qy), cosiq^qyll) — > cos{qnqJl) , cos{qn{qy + Qy)lf,) — > cos{qn{qx + Qx)lo)> cos{qnQylff) — > 
cos(g„Qa;;^), cos(q'„(Qy ± Q'y)lo) — > cos{qn{Qx ± Q'x)lo), 03(?n ± Qx'-'^) — > Qsiqn ± Qy'^'^), and cos{{qxQy - 
InQy^o) — > cos{{qyQx — qnQx)io)- Note that all other physical features of stripe phases (for example, the charge 
oscillation induced by the topological spin density of a skyrmion stripe phase in Appendix C and the perturbation 
theory for the stripe formation in Section IIIC, etc.) developed in the main text for the stripe aligned along y axis 
can also directly apply to the stripe phase along x axis in a similar way. The stripe phases constructed in a four level 
coherent trial wavefunction as shown in Section III B for even filling systems can also be obtained similarly. 



APPENDIX B: ANALYTICAL EXPRESSION OF THE FORM FUNCTION 



The explicit formula for the form function, ^ (q), defined in Eq. (7) can be evaluated analytically via a special 
function [17]. Here we just show the analytical results (for convenience, we define = (n(,,n^) and np = {n(j,n'^) 
to be the Landau level indices of a parabolic wide well system): 



nap.max^- ^'ap.max^- 



X exp 



cos^ 9{qyloY + {cos 9q Jo - sinOqJo)'^ Xl 



4Ai 
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sin'^ 9{qylof + (sin6'ga;^o + cosOqJpfXl 
AX2 



=p cos (^(ryy/cj j — /.{cusOqjJo — sin f^f/, /(, jAi 



V2xr 

cos^ Oiciylof + {cos 9q Jo - sm0qJo)'^Xl 
2X1 



V2X2 

,2 o/„ 1 \2 



sin^ 9{qylor + {sin9qJo + cos9qJo)^Xl 
2X2 



(Bl) 



where ± is the sign of n^^'' - n^^'^ for each bracket, n^^^^^^f^^^^^^''' = Min{Max) (n„^'\ n^^'^j , and m^^^'^ = 

\nj-''^ — np^'^\] Ai,2 = ^1,2/^0 = '^±/'^i,2 are dimensionless parameters, i™(a;) is the generalized Laguerre polynomial 
[61] and tan(26') = -2uj±_ijJ\\/ {ojI - o;^). 

As for a double well system with zero well width, we note that the noninteracting single electron wavefunction in 
Eq. (4) becomes Ly '^^'^ e'^^y'il:n\x + l'^k)y/5{z) for ojo +00, where ?/'i°-'(x) is the same as Eq. (5) with U replaced by 
the magnetic length, /q. Separating the z component and integrating q^ first, we obtain the following isotropic form 
function: 



A"" 



n. 



'aP,max 



exp 



g^j / ±qylo - iqJo 
' 4 J V ^/2 



2 



(B2) 



where q= \q±\, and all other notations are the same as above. 
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APPENDIX C: TOPOLOGICAL CHARGE DENSITY IN SKYRMION STRIPE PHASE 



According to Ref . [54] , the charge density induced by the topological isospin density in the double well system at 
u = 1 can be obtained from the isospin density function, ■m{f±): 

Ptopo{r±) = --e,j.,.rn{f±) ■ [di^m{fj_) x di,m{fj_)] , (CI) 

TT 

where the magnitude of m{rj_) has been normaUzcd to i. In this section we will show that the extra electron charge 
density, p^^ {r± ) , obtained in Eq. (42) for a skyrmion stripe phase is related to the induced charge density fluctuation 
calculated by Eq. (CI). Using the isospin density operator defined in coordinate space, Eq. (98), and the trial 
wavefunction in Eq. (36), we obtain the mean values of the isospinor in the 2D well plane after renormalization by 
the electron average density, {2ttIq)~^: 

{Mr±)) = ^ |l - 02(gx) [Aniqx, O, O) cos(g^(a; - Qyll/2)) + ^^^(fe, 0, 0) cos(g^(a; + Qyll/2))] | , (C2) 

and 

(Xf (rl)) = (J5(r-1))* =J2e3{q. - Qx)Am(9x, Q,, 0)e^(««-+«''^-^). (C3) 

Since Eq. (CI) is valid only to the leading order in (weak) isospin density modulation [54], we should take the limits 
of small isospin spiral and stripe orders (i.e. qlo <C 1 and qQylQ <C 1 for q = q^) in above equations and obtain 

{Iz{r±)) ~ ^ {cos(V'o) - 4A sin('^o) cos{qQyll/2) cos{qx)} = ^ cos{tp^/i2), (C4) 

where we have used tpk = '4'o + 4Acos(fcglQ) (Eq. (50)) for the last equation, and 

(J+(r-l)) = (J-(r±))* ~ ^^M(Q^,0)sin(V./,.)e^«^----'^ (C5) 

Note that the long wavelength limit behavior of (X±(rj_)) is different in wide and double well systems: for the former 

case, the two isospin states, and -H, have different Landau level indices, Hi and 7x2, and therefore A.ffij,{Q±,0) = 
AW^^{Qj_,0) ^ as \Qj_\ (see Eq. (Bl)). However, for the latter case, the two levels of different isospins are 

of the same Landau level index, A'^, and hence Afij.{Q±,0) = A^j^{Q±,0) ^ 1 as \Q±\ ^0. As a result, Eq. (C5) 
implies that the total magnitude of the isospin density in a wide well system is not a constant of winding vector, while 
it is a constant for a double well system. Therefore, in order to apply Eq. (CI) to the wide well system, we have 
to project the isospin vector onto the unit sphere by renormalizing its magnitude to 1/2 with a space independent 
constant and redefining the following isospin density vector, m(rj_): 

inx{r±) = i sin(^^/(2) cos((5j,a; + Qyy - 7), 
my{f±) = ^ sin(^^/;2) sin(Q^a; + Qyy - 7), 

m^(fj_) = i cos(?/'^./;2), (C6) 

which leads to the following charge density modulation via Eq. (CI): 

2 

Ptopo{r±) = --m{fj_) ■ dxm{fj_) x dym{r±) 

TT 



Qy ■ / I N dipk I 



47rlg ^-r-^/'o' Qk jk=x/ll 

r-^ QyAqsm{ilJo)sm{qx)/ir. (C7) 
In the same weak isospin density modulation limit, the extra local charge density shown in Eq. (42) becomes 
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Pex{r±) - [A^fr(gx,0,0)cos(9x(a;-(33;Zo/2)) -^M(9x,0,0)cos(gx(a; + QyZo/2))] 

~ -^Asin(V'o) [cos(g(a; - Qyll/2)) - cos{q{x + Qyll/2))] 

-QyAgsin(V'o)sin(ga;)/7r, (C8) 

which is the same as the charge density shown in Eq. (C8) above (but with an opposite sign). Such long wavelength 
charge density modulation is nonzero only when both the spiral order and stripe order are finite, and when the spiral 
winding vector is perpendicular to the stripe normal vector. Therefore this charge density modulation induced by the 
topological isospin density can be realized as a feature of isospin skyrmion stripe phase. 



APPENDIX D: RELATIONSHIP BETWEEN THE PERTURBATION STUDY OF THE STRIPE 
FORMATION AND THE COLLECTIVE MODE SOFTENING AT FINITE WAVEVECTOR 



In Section III C we discuss the suitability of the perturbation method to study the stripe formation if it results 
from an instability of a uniform phase via a two-step phase transition. In this section we will show that this method 
is equivalent to studying the stripe formation via the mode softening of a finite wavevcctor collective mode inside an 
isospin coherent or spiral phase, which is known as a standard method to investigate the stripe phase instability in 
a double well system [27,24]. Before we investigate the relationship between these two methods, it is instructive to 
mention that in general the collective mode dispersion can be obtained from the HF energy calculated by using a trial 
wavefunction similar to Eq. (36) but based on the symmetry-broken phase. More precisely, we can start from the 



following trial wavefunction, which is based on the rotated isospin basis, c 



1(2), fe 



obtained in Eq. (36): 
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Cl,fc 













gifeQxi^/2 cos(V;o/2) e-*'='3x/^/2 sin(V;o/2) 
_eikQ.ll/2 sin(V;o/2) e-*'=^-'o/2 cos(V^o/2) 









l,k-Qy/2 











(Dl) 



where we use Ci(^2) fc; V'Oj a-nd to denote the new state and the relevant variational parameters for this second 
rotation. For simplicity, we have taken the phase 7 = 7 = in the following discussion. As discussed before, we 
just need to consider a uniform ^0 to study the collective mode dispersion, i.e. when we take V'o to be small, the 
coefficient in front of leading order (quadratic of tpo) term of the HF variational energy calculated by the new many- 
body wavefunction, \^g), can give the dispersion of a low-lying collective mode. This method is basically equivalent 
to the conventional time-dependent Hartree Fock theory used in the integer quantum Hall system [17,33]. 

The study of stripe formation from the collective mode softening at finite wavevector inside the coherent phase region 
can be realized in the following isospin rotation picture: When the system has no coherence (i.e. the conventional 
uniform QH state), the isospin is polarized along isospin z axis. The trial wavefunction of a many-body phase shown 
in Eq. (36) can be understood as a result of a rotation in isospin space, tilting the isospin polarization from +z to 
+z (see Fig. 3(a)) with a tilting angle, tpo- If the spiral winding wavevector is not zero (i.e. a spiral phase), the 
isospin polarization vector changes its x and y components in the old isospin space with respect to the guiding center 
coordinate, k, but keeps its z component, (Iz), as a constant along the original z axis, i.e. no stripe order. When 
we apply a similar rotation, Eq. (Dl), onto the existing coherent or spiral phase, the isospin polarization vector 
prefers to have a second tilting angle, ipo, about the z axis (see Fig. 3(b)). As shown in Fig. 3(b), if the isospin 

vector winds about the new axis, z, (i.e. Qj^ ^ 0) the isospin projection onto the original z axis becomes periodically 
oscillating, showing a stripe structure in the original isospin space. The amplitude of such stripe phase is proportional 
to the amplitude of the second tilting angle, ipoj if it is small. Combining with the fact we mentioned above that the 
collective mode can be obtained from the HF variational energy by taking small Vo in the coherent phase, we can 
conclude that the stripe formation, if it results from a two-step second order phase transition, can be investigated by 
the mode softening of the collective mode at a finite wavevector inside the coherent (or spiral) phase. On the other 
hand, if this collective mode is always gaped at finite wavevector, we can conclude that there is no second order phase 
transition for a stripe formation, and the only possible stripe phase is that from the first order transition. 

To understand better the relationship between the method of studying the stripe formation from the collective 
mode softening at finite wavevector and the method of the perturbation theory developed in Section IIIC, we can 
apply Eq. (Dl) to Eq. (36) (assuming the first tilting angle, Vo Oi for the first rotation is nonzero and known by 
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minimizing the HF energy without stripe order in i/jq), and let the second tilting angle 'ijjo to be small. Using c^^^-^ ^ as 
an intermediate isospin basis, and considering the conventional Landau gauge, A[y]{f), with the second spiral winding 
wavevector being along x, i.e. Qj^ = (Qx,0), we obtain 



cos(^o72) - sin(^o72)(^o/2)e-''='3»'o) ^ 



Qy/2 



+ e 



-ikQ: 
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Similarly we can obtain the following energetically degenerate state by using Qj^ = {—Qx,0)' 
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Then we can define a new operator by adding above two equations together in the leading order of V'o: 

=1 L-»/cQx;^/2~t pikQ^ll/2-~"\ 

= eikQ.il/2 (cos(^o72) - sin(Vo72)(V'o/2)cos(A;Q,/2)j '^ik-Qy/2 

^^-ikQ.il/2 (sin(V;*/2) + cos(Vo72)(V!o/2)cos(fcQ,Zo')) 4,fe+Q,/2 
~ e^'^«'o/2 cos(Vfe/2)ct,,_Q^/2 + e-^'=««'o/2 sin(Vfe/2)4 .^q^/^, 



where 



V'fc = -00 + ^0 cos(fcgZo) 



(D3) 



(D4) 



(D5) 



for V'o ^ 1- Wc note that Eq. (D5) is exactly the same as Eq. (50) shown in Section IIIC (4 A = V^o). Therefore 
we have shown that studying the stripe formation instability by using Eq. (50) is equivalent to studying the finite 
wavevector mode softening of the collective mode inside a coherent or a spiral phase. Note that another kind of stripe 
phase with modulation along y direction can be obtained by the similar method using the other gauge, A[j,](r ), in 
which particle momentum is conserved along x direction. Therefore the perturbative method for the stripe formation 
developed in Section III C is justified and can be applied in a more general situation. 



APPENDIX E: A GENERAL STRIPE PHASE FUNCTION, ipK 

In general, we can use the following periodic function for tpk in Eqs. (36), (46) and (49): 

r Vi forO< |A:/g/a-m| <^(l-77)/2 

^^ = } -§;{k- ^/2) + Vo for ^(1 - r?)/2 < \kll/a - m\ < ^(1 + rj)/2 , (El) 
[ ip2 iov ^{1 + 7]) /2 <\kl^/a-m\ < 1/2 

where ipo = {4>i + 4>2)/2 and tp = ipi — ■^2- < ^, rj < 1, < ipi, V'2 < tt and m is an integer. The meaning of these 
variational parameters can be understood more clearly from Fig. 6. Note that when 77 = 0, we have a rectangular 
function, while for rj ^ 0, we have a smooth i/j^ to the linear oder. To calculate Qi{qn) {i = 1,2,3) for such ip]^, we 
can first evaluate the following quantity (let k' = kl^/a, q'^ = QxCl): 

Q^^Y.'^OsikqxlD^^'' 
^^'t' k>0 

[<(T--n)/2 ^{(l+»7)/2 ^ /■1/2 

= / dk'cos{k'q'^)e''^' + dfc'cos(A;'g;)e-*^('= -«/2)/€'7+iV'o + / cos(A;'g;)e^'^^ 

Jo Ji{l-n)/2 Ji{l+Ti)/2 
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/ /, \ • . / , XX sin(£(l — r?)o' /2) , ., . . . /sinfo' /2) sin(£(l + 7?)g' /2) 

= (cos(V'i) +ism(V'i)) — ^^u-ixi ^ + (cos(V'2) +isin(V'2)) ' 

+ 

+ 



2C7?cos(g;g/2) cos(V>o) cos((/./2) sin(^r?g;/2) - y cos(gr?g;/2) sin(y>/2)] 
2C7?sin(g;^/2) sin(V>o) [-y cos(y>/2) sin(e7?gy2) + ^m', cos(^r?g;/2) sin(y)/2)] 



. 24)?cos(g;C/2) sin(V;o) [-gr;g; cos(v?/2) sin(^r;g;/2) + ^ cos{im'J sin(y/2)] 
^.2^r?sin(gi^/2)cos(Vo) [-y cos( (^/2) sin(er?gy 2) + ^m'. oos{^VQ'J2) sin(y/2)] 

Then we have 



■'2„2„/ 2 2 ■ ^ ' 



2 1 1 

Q2{qx) = jrY. cos(fcgx/o) sin2(V'fe/2) = 1^ cos(fcfeZ^)(l - cos(V'fe)) = -5g,,o - iie[e], (E3) 

fe>0 fe>0 

and 

03(gx) = TT m cos(fcga;/o) sin('!/'/c/2) cos(V'fe/2) = XI cos(fcga;/o) sin(V'fc) = /m[e]. (E4) 
When replacing q'^ = q^a = 2-Kn, we have sin(g^/2)/(g^/2) = 5„^o- 



APPENDIX F: EVALUATION OF E 



i,o 



Here we show the analytical formula of the exchange energies used in Section VII for a double well system at 
y = AN + 2. For convenience, we let Qj^ = QiJo, P' = Plo, and d' = d/lo, q' = qlo to be dimensionless. We have 



-g 11 



2 r°° 



— 1 — p-^ c 

E'/{Q, 0) = ^ E yNN,NN{?i-) = -[^l dq' 

X i [l±Mq'P')e-<'''^' 
E^F°{Q±, 0) = 77- XI ^NN,NNi9±) cos{{Q^qy - Qyqx)ll) 



— J^ - (1 ± cos(g'(P^ cos e + P'y sin 6)) , 



= --/ dg'e-«V2 

'0 Jo 



-g /2 



-9 /2 



-'N 



-1 /'27r Ta 

2- y y (1 ± cos(g'(P^ cos + sin^)) e"?''^'] cos(g'(Q; sin0 - cos 9)) 



-0 ( 9 



/2' 



and 



±\ {jo {q'yJ{P^-Q'yy + {P; + Q'J') + Jo {q'yJ{P^ + Q'y? + {Pi-Q'.r)) 

-E^'°(0, an) = ^ X! ^iv'M-,JViv(9-i-) cos{qnqyll) 



(F2) 



^0 Jo 



-'N 



T J 2^1 YV^ cos(g'(P^ cos^ + P; sin^)) e-^'"^') cosCg^g' sin^) . (F3) 
Jo (g'^Pf + (P^ + gJ2) + Jo (g'^Pf + (P^-g;)2 



Jo(g'«;) ± 
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The most general expression of the exchange energy including both spiral and stripe orders {Ep°{Q±,qn) = 
1^1^ J2q^NN JViv(9J-) cos(g„Q'j,/o) cos(((5x% " Qy<lx)io)) Can be obtained similarly via the same strategy. 
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(intersubband) 
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n2 = (0,iV),s=T (i) 
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W2' 
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(intrasubband) 
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ni = (0,iV),s =i 


ni = (0,iV + l),s=t 


(skyrmion stripe) ^ 


(parity, spin & trans, 
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Dl 


double well 


4iV + l 


iV,Z = +l,s=t 


N,l = -l,s =T 


coherent 
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coherent stripe 
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parity & trans. 

in X (and y) direction 
parity & trans, in y 


D2 


double well 


4Ar + 2 




N,a = -l,s =T 


coherent^ 


parity & spin 



TABLE I. Table of the isospin notations and the many-body phases of different systems discussed in this paper, n = (n, n') 
is the orbital level index of the energy eigenstate of an electron in a parabolic well subject to an in-plane magnetic field (see 
Eq. (2)). s — ±1/2 and I = ±1 are spin and layer indices respectively, a = ±1 is for the symmetric(antisymmetric) state of 
the double well system. The many-body phases of Dl case are described in the layer index basis (sec Section VI), while they 
are described in the symmetry-antisymmetry (noninteracting eigenstate) basis in the D2 case. The in-plane magnetic field is 
fixed to be in +x direction. 

J: Such level crossing exists only when N > N* for N* = M.ax.{uJc/oj±,u)j_/oJc), see Ref. [43]. 

f: We do not really obtain a many-body phase within the HF approximation, but the experimental data of Ref. [15] suggests 
that a skyrmion stripe phase may exist in the W2' system (see Section V). 

§: = "commensurate canted antiferromagnetic phase" in the literature with spiral order in the layer index basis (see Section 
VII). 
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skyrmion stripe 


skyrmion stripe 
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spiral stripe 
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winding vector, Q±_ 
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(1) 
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isospin stripe 












normal vector, h 


-4— >■ 


(-) 


X 


I 


arbitrary 



TABLE II. A table summarizing the directions of the spiral and the stripe orders in different systems. The in-plane magnetic 
field is fixed in +x axis. For W2' system (wide well at even filling), we speculate that the skyrmion stripe is a possible ground 
state according to the experimetal results and our analysis in Section VII. Note that in a wide well system, the spiral order is 
degenerate at Ql = (0, ±Qy), while in a double well system, it is fixed along +y only in the commensurate state. Winding of 
the phases in Dl system is set by the in-plane magnetic field, and does not imply breaking of the translational symmetry (see 
discussion in Section VI D). 
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FIG. 1. Noninteracting Landau level energy spectra of a parabolic quantum well as a function of the parallel (in-plane) 
magnetic field, B\\. Wc choose following system parameters: B\_ = 3 T (ciJx = 5.2 meV), a;o = 7 meV, and \g\ = 0.44 for GaAs 
material, ni is the first Landau level index of n. Regions of intersubband and intrasubband level crossings are indicated by 
circles (see text). 
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FIG. 2. Schematic pictures of the noninteracting energy configuration in the regions of level crossing (or level near degeneracy) 
for the six systems discussed in this paper: (a) intersubband level crossing of a wide well bX v = 2N + 1, (b) intrasubband level 
near degeneracy of a wide well aX v = 2N + 1, (c) intersubband level crossing of a wide well aX v = 2N + 2, (d) intrasubband 
level crossing of a wide well aX v = 2N + 2, (e) interlayer level near degeneracy of a double well aX u = AN + 1, and (f) level 
crossing of a double well aX v = AN + 2. Solid(dashed) lines are for filled(empty) levels. The horizontal axis is the strength of 
in-plane magnetic field. Up(down) arrows denote the electron spin states. Note that the level indices for (a)-(d) are n = {n,n') 
for a parabolic wide well system, while they are (n, I) for (e) and (n, a) for (f) (see Table I) for different definitions of isospin 
indices. The circles in (b) and (e) denote the two near degenerate levels we consider. 



49 





(a) 



(b) 



FIG. 3. (a) Isospin polarization in a coherent or spiral phase: the new polarization axis is along z axis with an angle t/jo tilted 
from the original z axis. Finite winding wavevector causes {1x,y) oscillating in the real space but keeps (Xz) as a constant, (b) 
If the isospin coherent (or spiral) phase is unstable to form a new isospin spiral phase based on the rotated coordinate, the 
isospin winding about the new polarization axis, z, will cause a new periodically oscillating (T^) in the isospin projection onto 
the original z axis, showing a character of stripe phase. 
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FIG. 4. (a) A cartoon for the isospin spiral structure. Horizontal arrows denote the isospin direction and the spiral curve 
indicate the transverse isospin {Tx,Iy) order parameter, (b), (c) and (d) are the isospin coherent stripe phase, isospin spiral 
stripe phase, and isospin skyrmion stripe phase respectively. Shaded areas show isospin up (0) and down (®) domains, and 
arrows in the right hand side show the directions of By . Note that the isospin spiral direction is always perpendicular to _B|| . n 
is the normal vector of the stripes, which denotes the direction of isospin (Tz) modulation. For isospin coherent stripe in (b), 
there is no spiral order and therefore, in general, the stripe can have arbitrary direction with respect to the in-plane field. 
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FIG. 5. (a) A schematic noninteracting energy configuration for a general four level degeneracy of an even filling system. 
Level 1 is the highest filled level, level 2 is the lowest empty level, level 3 is the second highest level, and level 4 is the second 

lowest empty level as defined in Section IIIB. (a, s), in the right hand side denotes the quantum numbers of each level in the 
double well system {a is the parity quantum number and s is the electron spin), while (n^, s) is the quantum number for levels 
in the wide well systems (see Table I). Thick(thin) horizontal lines denote the filled(empty) levels. The thick upward arrow 
represents the density operator, p2,i, which annihilates one electron in level 1 and creates another one in level 2. This is the 
main mechanism and the order parameter for the many-body phases of the four level system near the level crossing region. The 
other upward arrows denote the density operators for pi, 3, p-i/s, Pa,?,, P4,i and P4/2 respectively from the left to the right, (b) 
The exchange interaction diagram for the coupling between p2,i and p4,3- The parity quantum number, a, has sign changed at 
the vertex. This is the main mechanism to stablize the canted antiferromagnetic phase in the double well system at i' = AN + 2. 




kln/a-m 



FIG. 6. A trial periodic function of the stripe phase function, Vfc, in Eq. (36). tpi, tp2, f, r/, and stripe period, a, are 
variational parameters (not all of them are independent). Its mathematical expression is shown in Appendix E. 
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FIG. 7. (a) Magnetoplasmon dispersion near intersubband level crossing region of a wide well system at v = 5 (system 
parameters are the same as in Fig. 1). The upper(lower) curves are for _B|| = 2.2 and 2.3 Tesla respectively. (The latter 

is calculated based on the isospin polarized basis and its negative energy at zero wavovcctor indicates a mode softening in 
a symmetry breaking phase. The correct curve for such plasmon mode should be modified based on the new coherent state 
as mentioned in the text.) Solid(dashed) lines are for dispersion along y and x axes respectively. The filled circle at g = 
denotes the energy of disconnected excitations, softening of which is a signature of parity symmetry breaking. Note that the 
energy of the circle is just slightly lower (< 0.01 meV) than the asymptotic magnetoplasmon energy in the long wavelength 
{q 0^). (b) Single particle energy as a function of in-plane magnetic field, of the same system. The dashed(dotted) 
lines are the spin up(down) levels of the two crossing isospin polarized states (including HF self-energy correction), while 
the solid lines are the energies of the many-body isospin coherent state, which breaks the parity symmetry of the system for 
2.2 Tesla < Bn < 2.40 Tesla. 
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FIG. 8. Magnetoplasmon dispersion of a wide well system at v — 1 with large in-plane magnetic field. The perpendicular 
magnetic field {B±) is 3 Tesla and the bare parabolic confinement potential (uio) is 3 meV. Solid, dotted, and dashed lines are 
for B\i = 20, 25, and 30 Tesla respectively. When By > 30 Tesla, the plasmon mode is softened at a finite wavevector in y 
direction (perpendicular to the in-plane field). The filled squares, triangles, and circles denote the energies of the disconnected 
excitation energy at gx = for By = 20, 25 and 30 Tesla respectively. 
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CO 0.2 0.4 0.6 0.8 

q (10*^ cm-') 

FIG. 9. Pcrturbativc energy, EpJ^t{q), for the stripe formation in a wide well system at ly = 1 in large magnetic field region. 
Dotted, dashed, and solid lines are for in-plane magnetic field, By = 31, 35, and 40 Tesla respectively. Thick and thin lines 
are obtained in Landau gauges, A[y]{f) and A[x]{f), respectively. This result clearly shows that a stripe phase can always be 
stablized to be along y direction (stripe modulation is in x direction) for By > 30 Tesla. System parameters are the same as 
used in Fig. 8. 
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FIG. 10. Phase diagram of a wide well system ai v = 1. System parameters are the same as used in Figs. 8 and 9 but with 
different confinement potential, wo, leading to different well widths. The well width is estimated from the size of the single 
electron wavefunction in the lowest subband of the parabolic confinement potential in the absence of in-plane magnetic field, 
i.e. 2(m*Ci;o)"^/^. 
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FIG. 11. Comparison of energies between three many-body phases of a wide well system aX v = Q near the level crossing 
point (Sji = 11.1 T) : spiral (solid line), spiral stripe (dotted line) and skyrmion stripe (dashed line). The system parameters 
are chosen to be the same as those in Fig. 1. The polarization is defined by 62(0) (see Eq. (37) and the text), which is zero 
for the isospin up state (spin unpolarizecl) and is one for the isospin down state (spin fully polarized). For all three many-body 
phases, the isospin winding wavevector Qx = —Q'±_ has been chosen to be 0.75 1^^ in y direction (perpendicular to the in-plane 
field) and zero in the x direction, which is the optimal value to minimize the HF energy if 62(0) 0, 1. (It can be also obtained 
from the wavevectors of the roton minimum in the magnetoplasmon excitations as calculated in Ref. [17].) The phases, 7 and 
7' in Eq. (49) are set to be zero. For the spiral stripe and skyrmion stripe phase, the stripe phase function, tpk, is calculated 
variationally by using equations in Appendix E and parameters defined in Fig. 6. For the convenience of comparison, here we 
have fixed V'o = 7r/2, ^ = 0.5, ?7 = and stripe period a = 0.67 ~ 10^'' cm, and only if is allowed to vary from — 7r/2 to +-k/2. 
When ip = Q, the energies of the stripe phases are the same as the spiral phase. (Note that the directions of the spiral stripe 
and the skyrmion stripe are different: the former one is along x direction while the latter one is along y direction, see Sections 
III A 2 and V.) We have checked that using other values docs not change the figure qualitatively, and also cannot stabilize any 
of these many-body phases. In our calculation shown in this figure, the lowest energy state is always the conventional QH state 
(i.e. either spin unpolarized or spin fully polarized states at 62(0) = or 1 respectively). However, we find that the energy 
differences between these states with the three many-body states (spiral, sprial stripe and skyrmion stripe) are very small, 
and the lowest energy of the stripe phase curve (in both spiral stripe and skyrmion stripe) is at finite ip, showing that if only 
a uniform spiral phase is stablizcd by some more sophisticated approximations, the stripe order may also be stablized with 
energy even lower than the uniform spiral phase. The stablization of a skyrmion stripe may be responsible to the resistance 
anisotropy observed in Ref. [15] in the strong in-plane field region (see text). 
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FIG. 12. Phase diagram of the bilayer system aX v = 5 m the presence of parallel magnetic field, (a) for fixed in-plane 
magnetic field (P/o = 27ri3||fi/o/'I>o = 1.5), and (b) for fixed tunneling amplitude (Asas = 0.013e^/eZo). In isospin language 
defined by layer index basis, phase I is an isospin coherent phase, phase II is an isospin spiral phase, phase III is an isospin 
stripe phase, phases IV and V are isospin spiral stripe phases. Note that phase V has a very long stripe period, and is related 
to the charge imbalance phase for the short-ranged interaction (see text). 
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FIG. 13. Phase diagram of double well systems aX v = AN + 2 in the presence of an in-plane magnetic field. The solid and 
dashed lines are for u = 2 [N = 0) and v = 6 {N = 1) without in-plane magnetic field, while the dotted lines are for u = 2 
with strong in-plane magnetic field {Plo = 15). The layer separation, d, is 0.067^0. FM is the ferromagnetic state with both 
filled levels being in the same spin polarized direction, and SYM is the symmetric state, where both spin indices are equally 
occupied in a symmetric orbital Landau level. In the middle is the canted antiferromagnetic phase, which breaks the spin and 
parity symmetries and is a commensurate state when finite in-plane magnetic field is applied. 




FIG. 14. Possible domain wall structure in integral quantum Hall system near level degeneracy region. Arrows denote the 
isospin polarization (see Fig. 3(a)). In the surface of each domain, isospin hybridization may be generated, resulting a coherent 
spiral order. The in-plane magnetic field is in +x direction. 
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